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Nuclear  Magnetic  Resonance  Saturation  and  Rotary  Saturation  in  Solids* 

by 

Alfred  G.  Redfield  f 

Division  of  Applied  Science,  Harvard  University 
Cambridge,  Massachusetts 

Abstract 

ZT  6>  3 

Nuclear  spin-lattice  relaxation  times  of  A1  in  pure  A1  and  Cu  in 
annealed  pure  Cu  have  been  measured  with  a nuclear  induction  spectrometer, 

by  the  method  of  saturation.  The  experimental  values  of  T,  are  4.  1 + .8 

Z7  & 3 ^ — 

milliseconds  for  Al  and  3.0  + .6  milliseconds  for  Cu  , in  reasonable 

agreement  with  theory. 

The  dispersion  mode  of  the  nuclear  resonance  was  also  observed,  and  it 
was  found  that%'(the  real  part  of  the  rf  susceptibility)  does  not  saturate  at 
the  same  level  as  the  absorption,  X",  but  remains  roughly  constant  out  to  a 
radio-frequency  field  intensity  of  about  2 gauss.  BothX'  and  %"  become 
narrower  and  nearly  Lorentzian  in  shape  above  saturation.  When  the  dc 

magnetic  field  modulation  is  increased  from  14  to  41  cycles  the  phase  of  the  ■ 

dispersion  signal  lags  behind  the  modulation,  presumably  because  the  modu- 
lation frequency  is  then  comparable  to  T,  . Large  dispersion  signals  above 

1 ■ 23 

saturation  have  also  been  observed  for  the  Na  resonance  in  NaCl. 

This  behavior  of  the  dispersion  mode  is  in  conflict  with  the  predictions 

of  Bloembergen,  Purcell,  and  Pound  and  of  the  Bloch  equations.  The  validity 

. 

of  these  theories  is  reexamined,  and  it  is  concluded  that  although  they  are 
applicable  to  nuclear  resonance  in  liquids  and  gases,  and  to  solids  at  small 
rf  intensities,  they  contain  incorrect  assumptions  as  applied  to  solids  at 
high  rf  power  levels  . The  theory  of  Bloembergen,  Purcell  and  Pound  is 
based  on  an  assumption  equivalent  to  that  of  a spin  temperature.  It  is  shown 
that  the  spin  state  cannot  be  strictly  described  by  a spin  temperature  because 
the  phases  of  the  spin  quantum  states  are  not  incoherent,  as  required  by  the 
temperature  concept.  The  transverse  decay  of  the  nuclear  magnetization 
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predicted  by  the  Bloch  equations  is  shown  to  be  partially  forbidden  by  energy 
and  entropy  considerations  if  a large  rf  field  at  the  resonance  frequency  is 
continuously  applied  to  the  solid. 

A theory  is  developed  which  is  applicable  only  to  solids  at  rf  magnetic 
field  intensities  well  above  the  saturation  level  and  which  is  in  reasonable 
agreement  with  the  experimental  observations.  The  Hamiltonian  is  trans- 
formed to  a coordinate  system  rotating  at  the  frequency  of  the  rf  field.  The 
resulting  time-dependent  parts  of  the  spin-spin  interaction  are  non-secular 
perturbations  on  the  time-independent  part  , and  can  therefore  be  ignored. 
Statistical  mechanics  is  applied  to  the  remaining  stationary  spin  Hamiltonian; 
specifically  it  is  assumed  that  the  spin  system  is  in  its  most  probable 
macrostate  (a  canonical  distribution  of  quantum  states)  with  respect  to  the 
transformed  spin  Hamiltonian.  This  assumption  is  justified  because  the 
transformed  spin  Hamiltonian  is  effectively  time  independent  and  the  spin- 
lattice  interaction  is  small,  and  it  is  analogous  to  assumptions  basic  to 
classical  acoustics  and  fluid  mechanics.  The  spin-lattice  interaction  merely 
determines  the  expectation  value  of  the  transformed  spin  Hamiltonian,  which 
can  be  readily  calculated  under  the  assumption  that  the  expectation  value  of 
the  spin  angular  momentum  of  each  spin  is  relaxed  independently  to  its 
thermal  equilibrium  value  by  the  lattice  in  time  Tj.  Both  fast  and  slow 
modulation  of  the  dc  magnetic  field  can  be  treated. 

"Rotary  saturation"  is  observed  by  applying  an  audio— frequency  magnetic 
field  to  the  sample  in  the  dc  field  direction  while  observing  the  dispersion 
derivative  at  resonance  with  a large  rf  field  H^.  When  the  audio-frequency 
approaches  -yHj  the  dispersion  signal  decreases  and  goes  through  a minimum. 
The  effect  is  easily  treated  theoretically  in  solids,  liquids  and  gas-es  by  using 
a rotating  coordinate  system,  and  is  a rotary  analogue  of  ordinary  saturation. 
It  is  a convenient  method  for  calibrating  rf  magnetic  fields  and  appears 
potentially  capable  of  providing  useful  information  on  the  solid  state. 
Experimental  data  on  rotary  saturation  are  presented  and  discussed. 

I.  Introduction 

This  paper  reports  an  experimental  and  theoretical  study  of  nuclear 
magnetic  resonance  in  solids  at  high  rf  magnetic  field  intensity.  Metallic 
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copper  and  aluminum  were  experimentally  investigated,  and  the  original 

objective  of  this  work  was  to  obtain  nuclear  spin-lattice  relaxation  times 

in  these  metals  for  comparison  with  the  observed*  Knight  shifts^  and  the 

3 

theory  of  Korringa,  which  relates  the  relaxation  times  and  Knight  shifts 
to  the  electronic  structure  of  the  metals.  The  spin  lattice  relaxation  times 
were  measured  by  the  method  of  saturation. 

In  the  course  of  these  measurements  it  was  found  that  the  dispersion 

mode  of  the  nuclear  resonance  signal  behaves  in  a way  which  is  in  conflict 

4 5 6 

with  the  existing  theories  ’ ’ of  magnetic  resonance  saturation.  As  a 
result,  the  validity  of  these  theories  as  applied  to  solids  was  reexamined, 
and  a theory  was  developed  along  somewhat  different  lines  which  appears 
to  agree  with  experiment  for  rf  magnetic  field  intensities  well  above  the 
saturation  level  This  theory  also  suggested  the  possibility  of  observing 
an  effect  which  we  call  "rotary  saturation,"  analogous  to  ordinary  saturation 
but  taking  place  in  the  effective  field  of  a rotating  coordinate  system. 

II.  Saturation  Data  and  Discussion 

4 

The  experimental  apparatus  was  a nuclear  induction  spectrometer 

7 8 

similar  to  those  previously  built  by  Weaver  and  by  Jeffries.  The  details 
will  be  described  elsewhere.  The  most  important  new  feature  of  this 
equipment  was  a suitable  calibrating  circuit,  permitting  relative  measure- 
ments of  rf  nuclear  susceptibility  to  be  made  independent  of  receiver  gain, 

rf  level,  and  other  variables.  The  output  of  this  spectrometer  yields  the 

9 

derivatives  of  the  real  and  imaginary  susceptibilities  X,'  and%."  The  rf 
intensity  was  determined  to  an  accuracy  of  better  than  5 per  cent  using  the 
method  of  rotary  saturation  described  below.  All  data  were  obtained  at 
room  temperature. 

Powdered  samples  of  pure  aluminum  and  annealed  pure  copper  were 
prepared  as  described  by  Bloembergen  and  Rowland.  ***  In  these  samples, 
electric  =quadrupole  effects  are  relatively  small  and  can  evidently  be 
neglected.  * 

The  relative  absorption  at  resonance  was  measured  by  integrating  the 
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recorded  absorption  signal.  The  absorption  data  are  plotted  in  Fig.  1, 

for  an  rf  frequency  of  7 . 6 Mc/s  and  a magnetic  field  modulation  frequency 

of  14  cps  . In  both  aluminum  and  copper  the  absorption  appears  to  follow 
4 5 6 

the  expected  ’ ’ dependence  on  rf  field  intensity  for  a system  of  dipolar 
coupled  spins: 

X"(vo,  Hj)  = X"(v0’  0)  [ 1 +<1/2)Y2  H*  Tx  g(v0>]  _1  , (1) 

where  'X."(v0>  Hj)  is  the  imaginary  part  of  the  nuclear  magnetic  susceptibility 

g(v)  is  the  shape  function  of  the  unsaturated  resonance,  v is  the  resonance 

°13 

frequency,  and  is  the  magnitude  of  the  rotating  rf  field. 

The  shape  function  g(v)  is  normalized  with  respect  to  integration  over 
frequency: 


g(v)  dv 


1 . 


(2) 


o 

Values  of  g(vQ)  were  obtained  by  integrating  the  integrated  absorption  signal, 
and  Tj^  was  obtained  from  the  data  of  Fig.  1 using  Eq . (1).  As  will  be 
discussed  below,  there  is  some  question  about  the  correctness  of  (1),  but  it 
is  believed  that  values  of  T j obtained  in  this  way  are  reasonably  accurate. 


Korringa  has  developed  the  theory  of  nuclear  magnetic  relaxation  and 
the  Knight  shift  in  metals.  His  theory  has  been  discussed  by  Bloembergen 
and  Rowland^  and  by  Holcomb  and  Norberg.  ^ Korringa  obtains  the  result 


,-l  4-rr  kT  +-.  r 87rV°gP  PF  ,2 

i = ~ h P (£o}  P(€o>f  3 1 


(3) 


where  vq  is  the  atomic  volume,  are  the  densities  of  electronic  states 

per  unit  volume  per  unit  energy  range  for  spins  up  and  down  respectively, 
evaluated  at  the  Fermi  level,  P^,  is  the  electronic  probability  density 
evaluated  at  the  nucleus  for  an  electron  in  the  Fermi  surface  (assuming 
that  the  electronic  wave  function  is  normalized  to  unity  over  the  atomic 
volume  vq) , g is  the  nuclear  g-factor  , andp  is  the  Bohr  magneton. 

In  deriving  this  equation  the  effect  of  electronic  correlations  and  of  possible 
P-character  (asymmetry)  in  the  electronic  wave  function  has  been  neglected. 


10 


.3  .4  .6 


.03  .04  .06  .08  .1  .15  .2 

H,  GAUSS 


Fig.  1.  Relative  absorption  at  resonance  X "(v0)  as  a 
of  rf  field  intensity.  The  units  of  X"(v  ) are  not  the  same 
aluminum  and  copper  in  this  figure. 


function 

for 
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The  latter  assumption  is  approximately  justified  unless  the  P-character 
part  of  the  electronic  wave  function  is  much  greater  than  the  S-character 
part . 

Korringa  also  obtained  a value  for  the  Knight  shift  which  can  be  written 
(AH/Ho)2  = h/ir  kTg2T1  , (4) 


where  is  given  by  (3),  AH  is  the  Knight  shift  in  gauss,  and  HQ  is  the 
applied  magnetic  field  in  gauss.  The  effect  of  electronic  correlations  is 
also  neglected  in  this  expression. 


In  Table  I are  listed  experimental  values  of  g(vQ),  T^  and  AH/Hq,  together 

with  the  theoretical  values  of  T ^ obtained  from  Eq . (4)  and  the  experimental 

values  of  AH/H  . Also  included  in  Table  I are  the  values  of  v P„m*/m 
o o F 

obtained  from  Eq . (3)  and  the  experimental  T^  values,  assuming  the  electrons 
in  the  metal  can  be  treated  as  a Fermi  gas  with  m*  = m,,  and  with  one  electron 
per  atom.  The  quantity  voPp.is  the  magnitude  of  the  square  of  the  wave  function 
at  the  nucleus  for  an  average  electron  at  the  Fermi  level,  relative  to  that 
predicted  by  the  plane-wave  approximation.  Since  relaxation  effects  due  to 
possible  P-character  in  the  electronic  wave  function  have  been  neglected  in 
(3),  the  values  of  vqPj,  m*  /m  given  in  Table  I must  be  regarded  as  upper  limits. 

Table  I 


g(vQ)  sec 


Tj  sec 


(b) 


AH/H 


(c) 


Tj  sec  (Theoretical) 

(m*/m)  v Pa. (b) 
o r 


Aluminum 
2 12  + . 2 x 10'4(a) 

4.1  + . 8 x 10"3 

. 237  x 10"  2 

5.06  x 10~3 

250 


Copper 

2 . 28  + . 2 x 10'4(b) 
3.0  + . 6 x 10'3 

. 161  x 10-2 
2. 26  x 10  3 
260 


a.  T.  J.  Rowland  (private  communciation) . 

b.  Present  research, 
c . Reference  1 . 
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The  discrepancy  between  the  theoretical  and  experimental  values  of 
for  copper  may  be  due  to  electronic  correlation  effects.  Similar  dis- 

12 

crepancies  have  been  observed  and  discussed  by  Holcomb  and  Norberg 

in  the  alkali  metals.  The  opposite  discrepancy  in  aluminum  may  be  due  to 

the  effect  of  possible  P-character  of  the  electronic  wave  function,  the 

existence  of  which  is  likely  because  the  valence  electron  of  an  aluminum 

atom  is  a P- electron.  In  copper  the  contribution  to  of  the  P-character 

jof  the  electronic  wave  function  is  probably  negligible.  The  data  in  both 

aluminum  and  copper  are  in  poor  agreement  with  extrapolations  from 
14  15  16 

measurements  ’ ’ of  T ^ at  low  temperatures. 

The  derivatives  of  the  absorption  and  dispersion  obtained  in  aluminum 
at  various  rf  intensities  are  shown  in  Fig.  2.  The  resonance  curves  observed 
in  copper  are  similar  at  high  rf  levels,  but  more  nearly  gaussian  at  low  rf 
levels . 

The  important  feature  of  these  curves  is  the  fact  that  the  dispersion 
does  not  saturate  at  the  same  level  as  the  absorption,  but  remains  roujhly 
constant  out  to  a level  of  several  gauss  rf.  As  will  be  discussed  below, 
this  behavior  is  in  strong  conflict  with  the  existing  theories  of  magnetic 
resonance  saturation  in  solids,  which  predict  that  both  the  absorption  and 
dispersion  should  decrease  (near  resonance)  at  the  same  level,  according 
to  (1).  Portis^  has  observed  similar  nonsaturation  of  the  dispersion  in  the 
case  of  paramagnetic  resonance  in  F-centers,  but  the  saturation  behavior 
of  the  absorption  was  different.  The  theory  developed  by  Portis  for  the  case 
of  F-centers  is  not  applicable  here,  because  the  different  nuclei  are  tightly 
coupled  by  the  dipolar  interaction. 

The  absorption  line  becomes  narrower  above  saturation,  as  previously 

17 

reported  by  Abell  and  Kmght,  and  similar  narrowing  occurs  in  the 
derivative  of  the  dispersion.  This  is  also  in  conflict  with  theory,  which 
predicts  that  the  absorption  and  dispersion  curves  should  both  broaden 
upon  saturation. 

The  audio  phase  of  the  signal  at  the  output  of  the  receiver  was  not  the 
same  as  the  phase  of  the  modulation  applied  to  the  dc  magnetic  field  Hq  . 

In  the  limit  of  small  modulation  frequency,  the  nuclear  resonance  signal 
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and  the  modulation  appear  to  be  in  phase,  and  the  output  of  the  lock-in 
amplifier  probably  represents  the  true  derivative  of  the  absorption  and 
dispersion.  As  the  modulation  frequency  is  increased,  however,  the  nuclear 
resonance  signal  lags  behind  the  modulation.  The  apparent  derivative  of  the 
absorption  or  dispersion,  as  indicated  by  the  lock-in  detector,  is  then  not  the 
true  derivative,  and  a resonance  signal  can  also  be  observed  by  changing 
the  lock-in  detector  phase  to  90  degrees  with  respect  to  its  normal  setting. 
This  phenomenon  is  evidently  not  directly  related  to  the  nonsaturation  of 
the  dispersion  mode  noted  in  the  paragraph  immediately  above.  The  phase 
lag  of  the  resonance  signal  relative  to  the  modulation  is  evidently  a conse- 
quence of  the  fact  that  the  modulation  period  is  comparable  to  Tj.  A 

similar  phase  lag  can  be  observed  in  systems  obeying  the  Bloch  equations, 

18 

as  experimentally  and  theoretically  investigated  by  Halbach. 

In  Figs.  3 and  4 are  shown  the  apparent  dispersion  derivatives  at 
resonance  in  copper  and  aluminum,  observed  with  two  different  modulation 
frequencies  and  with  the  lock-in  detector  adjusted  to  detect  signals  either 
in  phase  or  in  quadrature  with  the  modulation.  The  in-phase  data  obtained 
at  14^>are  apparently  very  nearly  what  would  be  observed  in  the  limit  of 
low  modulation  frequency,  ("slow  passage")  and  therefore  nearly  correctly 
represent  the  behavior  of  the  true  derivative  of  the  steady-state  dispersion. 


For  purposes  of  calibration  and  comparison,  Figs.  3 and  4 also  include 
the  derivative  of  the  dispersion  at  resonance  of  protons  in  water  at  the 
same  frequency  (thus  at  a smaller  magnetic  field)  as  that  used  to  obtain 
the  aluminum  and  copper  data.  The  water  was  heavily  doped  with  para- 
magnetic impurity  (manganese  sulfate)  and  the  relaxation  time  of  the  protons 
was  of  the  order  of  1/5  millisecond  . The  Bloch  equations  predict  that  well 
above  saturation  the  derivative  of  the  real  part  of  the  nuclear  susceptibility 
of  such  a sample  should  be  (at  resonance) 


BX' 

8H 

o 


X H T,„ 
o o 2H 

7T 

' 1H 


2Hj  T 


(6) 


where  % is  the  proton  contribution  to  the  static  nuclear  susceptibility  and 

° 4 

T and  T ^ are  the  transverse  and  longitudinal  relaxation  times  for  the 

protons  in  the  sample.  The  water  and  metal  samples  occupied  nearly  the 
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same  volume  in  the  nuclear  induction  head,  so  that  no  geometrical  correction 
is  necessary  for  comparison. 

The  quantity  plotted  in  Figs.  3 and  4 is  the  output  of  the  lock-in  detector 
divided  by  the  product  of  the  rf  input  level  to  the  nuclear  induction  head,  the 
modulation  amplitude,  the  total  number  of  resonated  nuclei  in  the  sample, 
the  Q-factor  of  the  receiver  coil  in  the  nuclear  induction  head,  and  the  over- 
all gain  of  the  receiver  and  lock-in  detector.  Of  these  corrections,  the 
first  and  the  last  two  were  given  automatically  by  the  calibration  circuit. 

In  the  limit  of  low  modulation  frequency  the  quantity  plotted  in  Figs.  3 and 
4 should  be  the  relative  contribution  per  nucleus  to  the  rf  susceptibility 
derivative  at  resonance.  The  accuracy  of  these  data  is  thought  to  be 
better  than  + 20  per  cent. 

Apparently  these  effects  are  not  limited  to  metals,  but  are  characteristic 

of  solids  in  general,  well  below  their  melting  points.  In  insulating  solids 

the  observed  dispersion  signals  are  in  most  cases  different  from  what  they 

are  in  metals,  because  the  relaxation  time  is  likely  to  be  longer  than  the 

period  of  the  dc  magnetic  field  modulation.  In  Fi&  5 are  shown  dispersion 

signals  obtained  from  two  NaCl  samples  of  different  purity.  In  the  impure 

NaCl,  the  relaxation  time  is  thought  to  be  long  compared  to  the  modulation 

period  and  short  compared  to  the  several  minutes  required  to  pass  through 

resonance.  A similar  dispersion  curve  can  be  obtained  from  distilled 

1 8 

water  at  high  rf  levels,  as  discussed  by  Halbach.  In  the  Harshaw  NaCl 
the  relaxation  time  is  clearly  of  the  order  of  the  time  (about  one  minute) 
taken  to  sweep  through  the  resonance,  as  indicated  by  the  asymmetry  and 
hysteresis  in  the  observed  dispersion  signal. 

We  conclude  this  section  by  noting  that  the  shape  and  magnitude  of  the 
observed  dispersion  signals  in  both  metals  and  insulators  can  be  accounted 
for  by  the  assumption  that  in  the  limit  of  large  the  Bloch  equations  hold 
with  Tp  T?,  rather  than  = 1/2  g(vQ)  as  is  required  to  yield  the  correct 
line  width  below  saturation.  This  assumption  will  be  more  or  less  justified 
in  the  next  section,  and  the  resulting  predictions  and  comparison  with 
experiment  will  be  discussed  there  and  in  section  IV. 


r 


(siiNn  Aavanayv)  Noisaadsia  ±N3avddv 


Fig.  5.  Apparent  dispersion  derivatives  in  samples  of  NaCl  of  different 
purity.  The  Na23  nucleus  was  observed  at  a frequency  of  about  7.5  megacycles. 
The  modulation  frequency  was  41~in  phase  with  the  lock-in  detector. 


III.  Theory 


A.  Previous  Saturation  Theories 

4 5 6 

We  now  discuss  the  previous  theories  ’ ’ of  magnetic  resonance  satu- 

ration which  lead  to  incorrect  predictions  of  the  magnitude  of  the  dispersion 
at  high  rf  field  intensities  in  solids.  Throughout  this  paper  we  are  concerned 
with  solids  in  which  the  nuclei  can  be  regarded  as  fixed  in  their  lattice  positions, 
and  we  neglect  diffusion  and  lattice  vibrations.  This  is  justified  only  well  below 
the  melting  point. 

5 

The  theory  developed  by  BPP  is  the  simpler  of  the  two  approaches  to 
saturation  in  magnetic  resonance,  and  is  based  on  three  assumptions.  The 
first  of  these  is  that  the  effect  of  the  spin-lattice  interaction  is  to  relax  the 
spins  to  their  equilibrium  state  in  a time  T^.  This  is  evidently  correct,  at 
least  to  a good  approximation,  in  most  cases.  The  second  assumption  is 
that  the  spins  interact  strongly  with  one  another,  so  that  energy  absorbed 
at  one  frequency  of  the  dipolar  broadened  resonance  line  is  quickly  trans- 
ferred to  all  the  spins,  whether  or  not  they  are  in  a local  field  exactly  corre- 
sponding to  the  applied  rf  frequency.  This  assumption  has  been  discussed  by 
20  6 

Bloembergen  and  by  Portis,  and  appears  to  be  justified  in  this  case  because 
of  the  possibility  of  mutual  spin  flips  between  neighboring  nuclei- brought  about 
by  the  dipolar  interaction.  The  third  assumption  of  the  BPP  theory  is  that 
the  complex  rf  susceptibility  is  proportional  to  the  difference  in  population 
of  the  nuclear  spin  levels  and  is  not  affected  by  the  presence  of  the  rf  field 
except  in  so  far  as  the  spin  level  populations  are.  When  the  spin  system 
becomes  saturated  the  populations  of  the  spin  levels  become  more  nearly 
equal  and  it  follows  directly  from  this  assumption  that  both  the  real  and 
imaginary  parts  of  the  rf  susceptibility  saturate  in  the  same  way. 

19 

If  the  line-shape  is  assumed  to  be  Lorentzian  below  saturation,  the 
predictions  of  the  BPP  theory  are  in  agreement  with  those  of  the  Bloch 
equations  (under  steady-state  condtio.ns)  and,  apparently,  with  experiment 
in  the  case  of  most  gases  and  liquids.  As  discussed  below,  the  BPP  theory 
also  predicts  the  correct  asymptotes  for  “Xi"  in  solids  for  the  limits  — *•  0 

and  >co,  and,  therefore,  the  correct  qualitative  behavior  of  %'•'  over 

-9- 
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the  whole  range  of  H,  (except  in  the  special  case  of  inhomogeneity  broaden- 
ing  ).  The  BPP  prediction  of  the  variation  of  Mz  just  at  the  onset  of 
saturation  is  also  correct.  Evidently,  these  predictions  are  valid  in  solids 
because  they  are  based  only  on  the  first  two  assumptions  of  the  BPP  theory. 
The  predictions  concerning  the  dispersion  above  saturation,  which  are  in 
marked  disagreement  with  the  experimental  data  of  section  II,  are  a conse- 
quence of  the  third  assumption,  and  it  appears  certain  that  the  trouble  lies 
here . 

20 

Bloembergen  has  pointed  out  that  the  BPP  theory  must  be  incorrect 
when  a large  rf  field  is  applied  to  a solid  in  which  T ^ 7^  T^  because  in  the 
limit  as  — >-oo  the  line  width  is  predicted  to  be  yH^  J T rather  than 

yH^  as  expected  from  the  uncertainty  principle  and  the  transition  probability 
associated  with  the  rf  field. 

The  assumption  that  the  complex  susceptibility  is  proportional  to  the 

difference  in  population  of  the  adjacent  nuclear  spin  levels  is  equivalent  to 

the  assumption  that  the  spin  system  behaves  as  if  it  were  at  equilibrium  at 

a spin  temperature  higher  than  the  lattice  temperature,  corresponding  to 

the  actual  distribution  of  nuclear  spins  among  the  21+1  levels.  Evidently, 

the  assumption  of  an  equilibrium  distribution  of  spin  states,  i.e.,  a spin 

temperature,  is  not  justified  in  the  case  of  saturation,  where  the  spins  are 

subject  to  a large  time-dependent  secular  perturbation.  Although  the 

amplitudes  of  the  spin  states  may  be  described  by  a Boltzmann  distribution, 

21 

the  phases  of  the  quantum  states  are  not  random,  as  would  be  the  case  at 
thermal  equilibrium  and  as  is  required  in  most  quantum-mechanical  pertur- 
bation calculations.  This  can  be  seen  from  the  fact  that  M and  M , the 

x y 

transverse  components  of  nuclear  magnetization,  are  not  zero  and  cannot 
in  general  be  made  arbitrarily  small  compared  with  M if  the  system  is 
saturated.  A random  distribution  of  quantum  phases  would  lead  to  zero 
transverse  magnetization.  For  this  reason  it  is  surprising  that  the  simple 
assumption  of  a spin  temperature  predicts  the  same  behavior  for  a Lorentzian 
line  as  the  more  detailed  assumptions  of  the  Bloch  equations. 

4 

The  Bloch  equations  are  based  on  the  assumptions  that  the  interaction 
of  the  spins  with  the  lattice  and  with  each  other  can  be  considered  inde- 
pendently of  their  interaction  with  the  externally  applied  magnetic  field. 
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andthatthe  effect  of  spin-lattice  and  spin-spin  interactions  can  be  described 

by  simple  relaxation  terms . They  predict  a Lorentzian  line  shape  and  saturation 

behavior  similar  to  that  predicted  by  BPP,  apparently  in  good  agreement  with 

experiment  for  most  liquids  and  gases.  Detailed  theoretical  justification  of  the 

22  23 

Bloch  equations  has  been  given  by  Wangsness  andBloch  and  by  Kubo  and  Tomita, 
but  their  work  is  applicable  only  to  liquids  and  gases  in  which  the  correlation 
time  of  the  motion  of  the  nuclei  is  short  compared  to  the  Larmor  period.  These 
authors  neglect  the  effect  of  the  rf  field  on  the  relaxation  process,  and  Bloch  and 
Wangsness  also  consider  the  spins  as  independent  of  each  other,  which,  as  they 
recognize,  is  not  justified  in  the  case  of  solids. 

The  Bloch  equations  have  always  been  recognized  as  a crude  approximation 

for  solids  well  below  the  melting  point,  since  they  predict  a L.orentzian  line 

19 

shape  rather  than  the  nearly  gaussian  shape  observed  experimentally . Thus  it 
is  natural  to  suppose  that  the  spin- spin  terms  of  the  Bloch  equations  (those 
involving  are  the  source  of  the  present  difficulty.  This  can  be  seeninmore 
detail  by  considering  a specific  solution  of  these  equations. 

Suppose  that  a circularly  polarized  rf  field  of  magnitude  an»jl  frequency 

(the  resonance  frequency)  is  continuously  applied  to  the  solid,  and  suppose  that 
initially  the  nuclear  magnetization  ^ is  in  the  direction  of  the  rf  field  and  of 
magnitude  Mj  Assume  that  Tj>;>  T^  as  is  the  case  in  most  solids  at  low 
temperatures.  The  terms  in  the  Bloch  equations  involving  T^  can  then  be 
neglected  during  a time  comparable  to  T^,  and  the  prediction  of  the  Bloch  equa- 
tions is  that  will  remain  parallel  to  and  will  decay  exponentially  to  zero 
in  a time  T^. 

In  the  course  of  this  decay  work  is  done  by  the  spin  system  on  the 

external  magnetic  field;  this  energy  can  come  only  from  the  internal  (spin- 
spin)  energy  of  the  spin  system.  The  energy  cannot  come  from  the  lattice, 
because  we  have  implicitly  neglected  the  spin-lattice  interaction  by  neglecting 
the  T^  term  in  the  Bloch  equations.  Speaking  somewhat  classically  and  loosely 
we  can  say  that  the  initial  state  corresponds  to  n excess  spins  aligned  in  the 
rf  field  direction  and  the  final  state  corresponds  to  no  excess  spins 
aligned  in  any  specific  direction  and  s/Y"  excess  spins  aligned  preferentially 
in  the  direction  of  the  magnetic  fields  of  their  neighbors.  Conservation  of 


TR206 


- 12- 


energy  requires  that 

nHj  =J/b  H,  (6) 

where  6H  is  of  the  order  of  the  r.m.s.  magnetic  field  at  a nucleus  due  to 
its  neighbors,  or  approximately  the  half-width  of  the  resonance  line  in  gauss. 

The  decay  predicted  in  time  T ^ by  the  Bloch  equations  is  an  irreversible 
process  in  a thermally  isolated  spin  system,  and  the  entropy  of  the  system 
must  increase.  The  initial  entropy  of  the  spin  system  is  the  same  as  that 
for  a spin  system  with  n excess  spins  parallel  to  a fixed  magnetic  field: 

S.  = S - kn2/N.  (7) 

SQ  is  the  maximum  entropy  of  the  spin  system  corresponding  to  complete 
chaos,  k is  Boltzmann's  constant,  N is  the  total  nutnber  of  spins,  and  for 
simplicity  we  assume  that  the  spin  is  1/2. 

By  analogy  with  (7),  the  entropy  of  the  final  state  is  expected  to  be 
approximately  given  by 

SF  = So  - ls/r2/N.  (8) 

This  assumption  can  be  justified  by  a detailed  calculation  similar  to  that  of 
Appendix  A discussed  below. 

From  (6),  (7),  and  (8)  it  follows  that  the  spin  entropy  increases  during 

r**' 

the  decay  only  if  H^<  6H.  Thus  for  sufficiently  large  the  Bloch  equations 
are  evidently  incorrect,  because  they  predict  an  irreversible  process  in 
which  the  entropy  decreases. 

The  actual  final  value  of  M after  such  a transverse  decay  can  be  estimated 
by  maximizing  the  entropy,  subject  to  the  condition  that  energy  must  be 
conserved.  We  assume  that  the  entropy  is  Sq  - k(n  +JY  )/N  and  require 

that  nH,  +JV 6 H be  a constant.  If  /V is  initially  zero  the  entropy  is  a maxi- 

L V 22-1 
mum  if  n,  and  therefore  M,  is  finally  (1  + 6H  /H,)  times  its  initial  value,  and 

2 2 A 

the  external  energy-M  is  H^/6  H times  the  internal  spin-spin  energy. 

To  summarize,  for  H^>  6 H the  spin  system  is  unable  to  take  up  the 
entire  energy  of  the  nuclear  magnetization  with  respect  to  the  rf  magnetic 
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field,  and  the  transverse  decay  predicted  by  the  Bloch  equations  is  partially 
forbidden.  As  will,  be  discussed  in  the  next  section  such  a process  is 
further  inhibited  by  the  rapid  decrease  of  the  suitable  transition  probabilities 
as  the  rf  field  is  increased. 

These  considerations  suggest  that  in  the  limit  of  large  the  x-y 

components  of  the  Bloch  equations  should  be  modified.  The  modification 

24 

can  most  easily  be  presented  in  a coordinate  system  rotating  about  the 
z-axis  at  the  rf  frequency  w (not  necessarily  the  resonance  frequency)  with 
its  x-axis  in  the  direction  of  the  circularly  polarized  rf  field  Hj.  In  such 
a coordinate  system  the  modified  Bloch  equations  are 

M = v[  M x H ] - M /T,  (9a) 

xr  ,L*wr  ~IerJxr  xr  2e 

M=y[M  xH  ] -M  / T,  (9b) 

yr  ll»r  erJyr  yr  2 

M = y[M  x H ] - (M  - M )/T.  (9c) 

zr  ,L'~r  — erJzr  ' zr  o 1 

The  subscript  r denotes  a quantity  measured  with  respect  to  the  rotating 

23 

coordinate  system,  and  H is  the  effective  magnetic  field.  H is  given 

gt  er 

by 


H = (H  - — ) 1 + H.x  , 

/w  e r o y ' •***  j-  l — r 

where  andjr  are  unit  vectors  in  the  rotating  system. 


(10) 


In  (9)  T^e  is  of  the  order  of  T ^ and  is  a transverse  spin-lattice  relaxation 
time,  since,  as  discussed  at  length  above,  the  spin-spin  collisions  are  unable 
to  relax  the  magnetization  in  the  direction  of  a large  rf  field.  T^  is  still 
the  correct  relaxation  time  to  use  in  the  y-component  of  (9)  because  a decay 
in  the  y-direction  (rotating  frame  of  reference)  involves  no  change  in  energy. 
Solution  of  (9)  leads  to  the  prediction 

lXoT2T2e 


X'  = 


<“-“</T2T2e+  *XT1T2+  ‘ 


*"=  X'(T2/T2e)W0-co) 


-1 


(Ha) 


(Hb) 


The  Lorentzian  dispersion  predicted  by  (11a)  is  in  reasonable  agreement 
with  both  the  shape  and  order  of  magnitude  of  the  observed  dispersion  at  high 


r 

I 
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rf  fields.  Unfortunately  (11)  is  quantitatively  applicable  only  at  the 
highest  rf  fields  attainable.  To  obtain  quantitative  results  at  lower  rf 
fields  it  is  necessary  to  consider  the  problem  in  more  detail. 

To  conclude  our  discussion  of  the  existing  saturation  theories  we  may 

4 

mention  another  defect  in  the  Bloch  equations.  As  usually  written  they 

predict  that  the  effect  of  the  spin-lattice  interaction  is  to  relax  the  nuclear 

A 

magnetization  toward  its  static  equilibrium  value  X^H^z^in  time  Tj. 
Actually  it  is  more  reasonable  to  assume  that  the  magnetization  relaxes 
toward  the  value  X^^H^corresponding  to  the  instantaneous  applied  field  fj. 
The  theory  and  experiment  in  the  present  paper  will  be  restricted  to  the 
case  so  that  H^^z^and^H  very  nearly  coincide  in  magnitude  and 

direction,  and  this  consideration  is  unimportant. 

B.  The  Rotating  Coordinate  Representation 

If  we  neglect  the  interaction  with  the  lattice,  the  complete  Hamiltonian 
of  the  spin  system  is^ 

X-,  = 8P H0 ^ !jz  + SPZHjlcoBut)  S IJx 


,2  ^ u; 


+ ® P k>j  ^ r jk  wj  ik  jk  ^~jk  £*jk*wk) 


k>J 


A.,  I. 
jk~J 


jk 

I, 


k ’ 


(12) 


13 


where  Hq  is  the  dc  field,  2H^  is  the  applied  linearly  polarized  rf  field  of 

frequency  w,  and  the  spin  operators  I.  are  expressed  in  units  of  U.  All  nuclei 

are  assumed  to  have  the  same  g-factor  and  spin;  the  case  of  two  magnetic 

ingredients  will  be  considered  later.  The  last  term  is  included  for  genfer- 

ality;  evidence  for  the  existence  of  such  a mutual  nuclear  interaction  in 

Z 6 

solids  has  been  discussed  by  Bloembergen  and  Rowland  and  by  Ruderman 

27 

and  Kittel . The  analogous  interaction  between  nuclei  in  molecules  was 

28 

first  suggested  by  Ramsey  and  Purcell. 


The  spin  wave -function  vJj  obeys  Shroedinger ’ s equation: 

= >fs^ 

We  will  use  the  transformations 


(13) 


TR206 


* = Ry0^r  = RyORa,t+ 


where 


R . = exp(  -iwt  / I . 

ZUt  JJ 


R = exp( -i®  y I.  ) 

y©  r'  *-y  jy' 


@ = tan  ^[yHj/(w-<*)  Q)] 


Then  it  follows  that 


-ih  9ib  /8t  = (-TxjjY7!  . + R ,^£r  ~1)\b  = ^6  iL 

A-1  zi  zwt  S zwt  ,Yr  errr 


-ihSib  /3t  =Rve^  R e'1  +.  s >£„+ 


The  unitary  transformation  R-Z(jjt  can  be  regarded  as  a transformation 
to  a coordinate  system  rotating  about  the  z-axis  with  frequency  w (note  that 
this  is  the  rf  frequency,  not  necessarily  the  resonance  frequency  tj0o).  As 
above,  the  subscript  r is  used  to  denote  quantities  in  this  coordinate 
system  and  J^r  is  the  effective  Hamiltonian  in  the  r-system. 


Ry0  corresponds  to  a further  transformation  to  a coordinate  system  fixed 


with  respect  to  the  r-system,  whose  z-axis  coincides  with  the  effective 
magnetic  field  direction,  and  whose  y axis  coincides  with  the  y axis  of  the 


r-system.  0 is  the  angle  between  the  effective  magnetic  field  H in  the 

24 

rotating  coordinate  system  and  the  z-axis  in  the  fixed  system.  Quantities 


in  this  system  are  denoted  by  the  subscript  p. 


The  Hamiltonians  Hr  and^^p  in  (18)  and  (19)  are  readily  obtained  from 


the  fact  that  the  operators  I.  transform  like  vectors  under  the  rotations 


corresponding  to  R and  R _ . Since  the^I,.  occur  only  in  scalar  products 


the  same  result  can  be  achieved  by  applying  the  inverse  rotation  to  the 


other  vector  of  the  scalar  product.  Thus  I.*  I.  is  invariant  under  these 

/w  K 


transformations  while  the  r ..  are  transformed  under  the  inverse  rotation. 

jk 


This  procedure  is  simply  an  algebraic  shortcut;  actually  thejr^  are  in- 


variant parameters  and  the^  undergo  transformation  In  this  way  we  get 


I 


l 


1 


T r 
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- gPHer  • Zij  + X (Ajklj  -Jk+  Bj^,)  (20) 

j ^C>  J 


+ time -dependent  part, 


' «f>H.r  Z IJz  + .E  ' ik  + VeVks 

4 Djkpt  IJ+Ik++  'j-'k-1  + EjlvC<IJ+  + Ij->Ik*+IJ*(Ifc++  Ik->1} 


+ time -dependent  part, 


where 


Ajk  = Xjk+g2p2rIk(3/2  1/2> 

Bjk=-3gVr:k3(3/2  42.k-l/2) 


jk 


Ajkp=  Ajk  +(3/2co82®-l/2)(A.k  - Xjk) 


Bjkf  = (3/2co»^  O - 1/2)  Bjk 


Djkp=1/45i"  ®Bjk 


Ejj^  =1/2  sin  0 cos  © B^k 


H 


er 


ft* 


= j H2  + (Ho  - u/y  Y 


(21) 

(22) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 


Here  I.+  are  the  raising  and  lowering  operators  I.  + il.  having  selection 
J ^ J J/ 

rules  aMj.  = +1  respectively,  and  £,jk  is  the  z-direction  cosine  ofjTjk-  The  time 

dependent  parts  of  these  Hamiltonians  contain  terms  like  Ij  ^Ik_,  lj+Ikz*Ij+Ik+, 

etc .,  wjLth  time  dependences  exp(  + iwt),  exp(+2iwt).  + 

We  now  assume  Ho»Her=6H  and  apply  Dirac  perturbation  theory  to 
Schroedinger 1 s equation  in  the  rotating  coordinate  systems  (equations  18  and  19). 
The  time -dependent  part  of  ■^r(<:>r^!^p)  i®  a nonsecular  perturbation  on  the  time- 
independent  part.  It  connects  states  for  which  aMj  = + 1 for  one  or  two 
nuclei,  and  for  it  to  be  secular  these  states  must  differ  in  effective  energy 
by  + liw  - + ■y'h  Hq  or  + 2hw  = 2y^H0*  Actually  such  states 


29 


.r<2j 


J 
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will  differ  in  effective  energy  by  about  c-yti(Her  + 61^),  with  c =0,  1 , or 

2,  so  the  condition  for  a secular  perturbation  is  not  satisfied  and  the  time- 

30 

dependent  part  of  can  be  ignored  for  most  purposes.  It  will  be 

noted  that  these  time-dependent  terms  are  the  same  ones  neglected  by 
25 

Van  Vleck  in  his  calculation  of  the  moments  of  the  absorption  line,  plus 
a term  corresponding  to  the  component  of  applied  rf  field  rotating  in  the 
opposite  direction  from  that  of  the  rotating  coordinate  system. 


If  Hgr  is  very  large  (but  still  much  less  than  Hq)  the  fourth  and  fifth 

terms  of  7s/  n ar€f  nonsecular  perturbations,  since,  being  time-independent 
ep 


ep 


to  be  secular,  they  must  connect  eigenstates  of  the  first  three  terms  of 
having  the  same  effective  energy,  whereas  actually  they  connect  states 
differing  in  effective  energy  by  about  one  or  two  times  Y^(Her  + 6H).  Under 
these  circumstances  we  can  ignore  these  terms  and  %(-  conserves  M , 

° ep  zp 

the  nuclear  magnetization  in  the  effective  field  direction.  It  should  be  noted 
that  to  ignore  these  terms  the  transition  probabilities  associated  with  them 
must  be  much  smaller  than  T ^ . Thus  this  procedure  is  valid  only  for 

Her  very  much  larger  than  the  resonance  line  width. 


Conservation  of  M leads  to  a system  of  equations  similar  to  (9),  whose 

z p 

solution  in  the  limit  of  large  is  identical  with  (11).  We  do  not  go  into 
details  because  the  theory  is  not  applicable  to  any  experimentally  attainable 
situation.  Experimentally  we  are  limited  to  rf  field  intensities  comparable 
to  the  resonance  line  width. 


C.  The  Canonical  Distribution  in  the  Rotating  System 

In  the  previous  section  we  obtained  a transformed  spin  Hamiltonian 

which  was  effectively  independent  of  time.  A time-independent  Hamiltonian 

is  convenient  to  work  with  because  the  concepts  of  statistical  mechanics 

29 

can  be  more  easily  applied  to  it  and  because  the  effective  energy  of  the 
spin  system  can  change  only  through  the  spin-lattice  interaction.  The  spin- 
lattice  interaction  terms  of  the  Hamiltonian  acquire  an  additional  explicit 
time  dependence  as  a result  of  the  transformation  to  the  rotating  system, 
but  if  is  large  enough,  these  terms  can  be  considered  as  a small 
perturbation  which  transfers  effective  energy  between  the  spin  system  and 
the  lattice  and  eventually  determines  the  value  of  the  effective  energy  of  the 
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spin  system. 

As  in  the  statistical  mechanics  of  stationary  systems  we  consider  a 
Gibbsian  ensemble  of  systems,  each  consisting  of  the  solid  with  the  rf 
field  continously  applied  to  it.  Suppose  that  the  (time-dependent)  spin- 
lattice  interaction  is  somehow  turned  off,  and  that  the  average  expectation 
values  of  the  lattice  Hamiltonian  and  transformed  spin  Hamiltonian  ('Her) 
are  known.  In  the  absence  of  further  information  the  distribution  of  states 
in  the  ensemble  is  still  highly  ambiguous,  but  the  most  probable  distribution 
is  a canonical  distribution  of  states  with  respect  to  the  lattice  and  the  trans  - 
formed  Hamiltonian.  Associated  with  this  canonical  distribution  are  two 
temperatures,  determined  by  the  canonical  average  expectation  values  of 
the  lattice  and  transformed  spin  Hamiltonians . One  of  these  is  the  lattice 
temperature,  which  is  presumably  positive  in  any  physically  attainable 
solid.  The  other  will  be  called  the  effective  spin  temperature,  which  is 
positive  or  negative  depending  on  whether  the  average  expectation  value 
of  the  effective  spin  Hamiltonian  5^er  is  negative  or  positive.  The  situation 
is  closely  analogous  to  that  in  the  fixed  coordinate  system  if  the  spin-lattice 
interaction  is  turned  off  and  if  no  rf  field  is  applied.  In  that  case  the  most 
probable  state  of  the  system  is  also  represented  by  a canonical  distribution 
of  states,  with  lattice  and  spin  temperatures  which  are  not  necessarily  equal. 
Negative  spin  temperatures  are  easily  attained  (with  pulse  techniques)  in 
insulating  crystals  in  which  the  spin-lattice  interaction  is  very  small.  In 
both  the  fixed  and  the  rotating  coordinate  systems  the  lattice  and  (effective) 
spin  temperatures  can  be  different  because  the  spin  and  lattice  Hamiltonians 
commute,  and  it  is  sensible  to  talk  about  canonical  distributions  of  states 
only  because  the  lattice  and  (effective)  spin  Hamiltonians  are  time-independent 
(except  for  some  nonsecular  perturbations)  in  both  cases  . 

We  now  consider  what  happens  if  we  turn  on  the  spin-lattice  interaction 
and  wait  for  a time  long  compared  to  the  spin-lattice  relaxation  time,  but 
short  compared  to  the  time  required  for  the  rf  field  to  heat  up  the  lattice 
appreciably.  In  the  absence  of  the  rf  field  and  in  the  fixed  coordinate 
system,  the  spin  system  will  approach  a canonical  distribution  of  states 
with  equal  spin  and  lattice  temperatures.  In  the  presence  of  the  rf  field,  in 
the  rotating  system,  the  spin-lattice  interaction  will  change  the  average 


■ 
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expectation  value  of  the  effective  spin  Hamiltonian  to  some  quasi -equilibrium 
value.  We  have  no  rigorous  assurance  that  the  spin  system  will  remain  in 
a canonical  distribution  of  states  with  respect  to^fer  but  we  assume  that 
it  will.  This  assumption  can  be  regarded  as  an  admission  of  our  ignorance 
concerning  the  system;  lacking  detailed  information,  we  simply  assume 
that  the  system  is  in  its  most  probable  state  for  the  limit  of  zero  spin- 
lattice  interaction.  The  actual  effective  spin  temperature  in  the  steady 
state  is  determined  by  the  lattice  temperature  and  the  spin-lattice  inter- 
action, and  depends  on  Hq,  H^,  and  w.  Unlike  the  static  case  with  no  rf 
field,  the  effective  spin  temperature  is  different  from  the  lattice  temperature 
and  can  actually  be  negative.  The  reason  for  this  difference  is  that  the 
spin-lattice  interaction  in  the  rotating  system  contains  an  explicit  time 
dependence  which  is  not  present  in  the  fixed  system. 

An  analogous  problem  is  that  of  a gas  in  a fairly  well-insulated  bottle, 
connected  to  one  or  more  temperature  baths  by  heat  leaks.  To  find,  say, 
the  pressure  of  the  gas  it  is  necessary  to  assume  immediately  that  the 
gas  molecules  are  in  their  most  probable  state  (a  Boltzmann  distribution) 
subject  to  the  constraint  that  they  have  some  definite  energy  (corresponding 
to  the  gas  temperature) . The  problem  is  then  reduced  to  finding  the  gas 
temperature  as  determined  by  the  various  heat  leaks  (corresponding  to  the 
spin-lattice  interaction)  and  temperature  baths  (corresponding  to  the  lattice). 

If  the  bottle  is  constrained  to  move  it  is  necessary  to  transform  to  the 
bottle's  coordinate  system  before  applying  statistical  mechanics  to  the  gas, 
in  analogy  to  the  rotating  coordinate  transformation  used  here.  All  of 
classical  acoustics  and  fluid  mechanics  are  based  on  assumptions  similar 
to  those  used  in  this  paper  since  it  is  always  assumed  that  matter  possesses 
the  same  thermodynamic  properties  viewed  from  a suitably  moving  co- 
ordinated system  and  in  suitably  small  pieces  as  it  does  at  rest  in  a fixed 
coordinate  system.  When  temperature,  pressure,  or  velocity  gradients 
in  a gas  become  too  large  this  assumption  breaks  down  (i  . e , at  low  i 

pressures  and  in  shockwaves)  and  the  theory  becomes  difficult.  Like- 
wise in  the  case  of  spins  when  the  spin-lattice  interaction  becomes  too 
large  the  assumptions  used  here  break  down  and  the  theory  is  difficult 

A rigorous  justification  of  this  procedure  would  be  extremely  difficult, 

•i 

j 
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if  not  impossible.  It  might  be  based  on  something  analogous  to  the 
Boltzmann  H-theorem.  Thd  least  that  can  be  said  for  thia  assumption 
is  that  above  saturation  it  is  the  only  simple  one  that  is  not  obviously 
wrong . 

The  assumption  that  the  spin  system  is  in  a canonical  distribution  of 

states  with  respect  to  2feT  leads  immediately  to  the  conclusion  that  the 

magnetization  is  in  the  direction  of  the  effective  field  ^er»  just  as  in 

the  static  case  the  magnetization  M is  in  the  direction  of  The  effective 

external  energy  of  the  spin  system  is  • ^Her  =.■*•  ^zp^er5  where  i® 

the  z p component  of  magnetization  in  the  p coordiante  system.  We  define 

a spin-spin  energy  in  the  rotating  system  which  contains  all  the  spin-spin 

terms  of  : 

er 

Xss‘  j^L'Vj  -Jk  + Bjki.W  <29> 

The  internal  spin-spin  energy  is  comparable  to  the  external  energy  when  the 

effective  external  field  H becomes  comparable  to  the  local  fields  at  the 

er 

nuclei  due  to  their  neighbors  (i . e . , approximately  the  line  width  in  gauss) 

31 

just  as  in  the  static  case.  In  Appendix  A it  is  shown  that 

(30) 

where  denotes  the  canonical  average  expectation  value  of  the  operator 

0,  and 


6h2=4<ahz>  + yr 


TT 
AV  b 7 N 


k>j 


A2  . 


(31) 


Here  <AH  V v is  the  second  moment  of  the  unsaturated  resonance  line  as 
A V 25 

calculated  by  Van  Vleck. 


The  reasoning  which  leads  to  (30)  is  essentially  the  same  as  that  used  in 

Section  111(b)  to  find  the  state  of  maximum  entropy  of  the  spin  system  after  a 

2 

transverse  Bloch  decay,  and  the  6H  used  there  is  the  quantity  given  by  (31). 
The  state  of  maximum  entropy  of  a system  for  a given  energy  is,  of  course, 
described  by  the  canonical  distribution  of  states  corresponding  to  that  energy. 

The  problem  is  solved  if  the  value  of«^fer>^  can  be  determined,  because 
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we  can  write 


«J¥er»=  -MZpHer  + «#ss» 

* -MzpHert‘+  6H2/Hfr). 


(32) 


Transforming  back  to  the  fixed  system  we  have  (since  is  in  the  direction 
of  H ) 

/W  er 


M — M „ sin  ® cosut, 
x zp 

(33a) 

M ^ M sin  © S' n wt, 

y zp 

(33b) 

M si  M cos  0 - 

z z p 

(33c) 

To  determine  the  expectation  value  of  ^=£er  we 

use  a simple  relaxation 

assumption  to  account  for  the  spin-lattice  interaction.  The  physical  reason- 
ing which  follows  is  justified  in  more  detail  and  under  more  general  con- 
ditions in  Appendix  B.  We  assume  here  that  the  effect  of  the  spin-lattice 
interaction  is  to  relax  each  nucleus  independently  into  its  equilibrium  state 
in  a time  T 


[ a/at]SL«j,.»  = -«<i>>  - io)/ Tr 


(34) 


where  the  left-hand  side  is  the  spin-lattice  contribution  to  the  time  derivative 
of  the  expectation  value  ofj^.,  andjj^  *s  static  thermal  equilibrium  value 
ofJL,  given  by 

Jo  =j  gpl(I+l)H/kT.  (35) 


Here  the  applied  field  H can  be  closely  approximated  by  H 


A 

Z 

O 


( since 


It  might  seem  more  logical  to  assume  that  the  nuclei  are  relaxed  along 

the  effective  field  instead  of  the  actual  applied  field,  so  that  we  should  use 

H in  (35)  instead  of  H.  That  this  would  be  incorrect  follows  from  the  fact 
»er  -w 

that  the  electrons,  which  are  responsible  for  the  relaxation,  are  almost 
completely  unaffected  by  the  rf  field,  and  effectively  see  only  the  large  field 
Hq  in  the  z-direction.  Furthermore,  if  we  were  to  use  the  effective  field 
in  (35)  we  would  get  obviously  incorrect  predictions  for  the  dispersion. 

Equation  (34)  implies  that  T^  , the  transverse  electronic  relaxation 
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time  defined  previously,  is  equal  to  T,.  This  assumption  is  in  accord  with 
12  23  32  X 

theory  ’ ’ for  metals,  assuming  the  electrons  are  in  a Fermi  distri- 

bution of  states . 

We  assume  that  the  spin-lattice  interaction  does  not  perturb  the  spin 
system  cononical  distribution  appreciably  except  to  bring  about  a slow  change 
in<C<C^’er^)>  ^ » the  transformed  spin  Hamiltonian  average  expectation 

value.  The  spin  lattice  relaxation  can  be  fictionally  regarded  as  a two-step 
process.  The  first  process  is  the  scattering  of  the  nuclear  spin  into  a 
completely  random  orientation,  in  a time  T^.  The  second  process  is  a 
scattering  of  the  spin  into  an  orientation  with  probability  governed  by  the 
Boltzmann  distribution  of  states  with  respect  to  the  externally  applied 
magnetic  field,  in  an  infinitesimal  time  after  the  first  scattering.  These 
two  processes  correspond  respectively  to  the  two  terms  on  the  right-hand 
side  of  (34).  The  change  ira(<^£er^with  time  is  the  sum  of  the  separate 
changes  brought  about  by  these  two  processes. 

The  rate  of  change  of«^>fer»  due  to  the  first  (random)  scattering  is 

W-%  »=  + MzpH er/T,  - 2 «MSS»  /Tr  (36) 

The  external  effective  spin  energy  -M^H^  (expectation  value  of  the  first 
term  °^^er)  proportional  to  the  sum  of  expectation  values  of  components 
of  the  1^^,  and  is  thus  expected  to  decay  to  zero  in  time  T ^ for  random 
scattering  of  the  spins,  corresponding  to  the  first  term  in  (36).  The  spin- 
spin  energy<<J=£ _ on  the  other  hand,  is  quadratic  in  the  components  of 

DO 

the^Ij,  and  is  therefore  expected  to  decay  at  twice  the  relative  (logarithmic) 
rate  of  the  external  energy;  thus  the  factor  two  in  the  last  term  of  (36). 

The  effect  of  the  second  (Boltzmann)  scattering  is  to  change  the  nuclear 
magnetization  at  a rate  +M  /T,,  where  M = X H—  XH  £,(in  the  fixed  co- 
ordinate  system).  The  corresponding  rate  of  increase  of  the  effective 
external  energy  -M-H  is  -(M  /T,)-H  . The  change  in  the  spin-spin 

energy  is  negligible,  since  the  local  fields  at  the  nuclei  are  random  in 
orientation,  to  a very  good  approximation.  Thus  the  rate  of  change  of 

due  to  the  Boltzmann  scattering  is  the  same  as  that  of  the  effective 
external  energy: 


v 
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t8/8tlsLB«^Cr»  = 'MoHer  cos  ®/Tl 


(37) 


To  obtain  the  steady-state  value  of  Mz^we  use  (30)  and  (32),  and  set 
the  sum  of  (36)  and  (37)  equal  to  zero.  The  result  is 


M 


M cos  O 
o 


(38) 


op  1+26H*7H^ 

er 

where  is  the  quasi -equilibrium  value  of  for  constant  Hq,  H^,  and 

CO  • 

The  dispersion  is  given  by 


M 


X' 


xr 


2H, 


M sin  0 
OP 

2H, 


(39) 


or 


X'  = 


M0Y(“0-w) 

2[  (Wo-co)Z+Y2(Hj  + 26H<i)] 


(40) 


so  that  above  saturation  the  dispersion  is  Lorentzian  and  the  dispersion 
derivative  at  resonance  is 


ax  _ 

9Ho 


X H 
o o 


2(Hj  + 26H  ) 


(41) 


Equation  (33) implies  that  the  nuclear  magnetic  absorption  is  zero  above 

saturation.  Actually,  (33)  is  an  approximation  and  Mr  is  not  precisely  in 

the  direction  of  H . At  any  rf  level  there  is  finite  absorption  which  can 
~~er 

be  predicted  by  invoking  conservation  of  energy  in  the  fixed  system.  Energy 

is  transferred  from  the  spin  system  to  the  lattice  at  a rate  — Ho(Mz-Mq)/T  ^ • 

Equating  this  to  the  energy  absorbed  by  the  spin  system  from  the  rf  field, 

2 

which  is  we  Set 

X"  = Ho(2H1X'  cot  0 - M^)/2T1H^  u (42) 


} 
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or 


X"  = 


y2(H?  + 26H2)H 


2[(  cj  - w Q)' 


+ y“(H7+  26H 


M 

o o 

2)]t 


~TP~ 

1H1W 


(43) 


This  is  a Lorentzian  line,  of  the  si.me  width  as  would  be  predicted  by  the 

2 2 2 -1 

Bloch  equations  below  saturation  for  T ^ = \ (H^  + 26  H ) . At  resonance 

(43)  agrees  with  the  asymptotic  values  of^G'  given  by  the  Bloch  and  BPP 
theories.  The  same  is  true  of  off-resonance  in  the  limit  of  very  large  H^. 
This  agreement  corresponds  to  the  fact  that,  for  any  theory,  under  these 

conditions  M — *-0  and  %"  is  uniquely  determined  by  this  simple  conversion 

z 

of  energy  argument.  The  fact  that  all  theories  yield  the  same  asymptotic 
2C"  means  that  the  saturation  method  of  obtaining  T^  used  in  Section  II  is 
substantially  correct.  This  is  not  necessarily  true  if  the  maximum  ab- 
sorption derivative,  rather  than  the  integrated  absorption  derivative,  is 
used  to  obtain  T^.  Error  in  the  saturation  determination  of  T j can  also 
be  introduced  by  the  use  of  a magnetic  field  modulation  period  comparable 
with  or  less  than  T^,  as  is  usually  the  case  experimentally.  Fast  modulation 
is  treated  in  Section  M(e). 


The  theory  developed  above  is  valid  only  for  large  Hgr . If  Her  becomes 
comparable  to  the  value  of  at  which  the  absorption  begins  to  saturate, 
the  spin-lattice  interaction  can  no  longer  be  regarded  as  a small  pertur- 
bation and  the  magnetization  M will  tend  toward  the  z-axis  with  a value  M , 

,w\  O 

rather  than  toward  the  direction  of  H . Treatment  of  the  intermediate 

er 

case  would  be  extremely  difficult. 


The  theory  is  also  invalid  in  the  case  of  large  Hej.  but  small  H^;  i.e. , 

in  the  case  of  ®— 0 or  (off-resonance).  In  this  case  the  spin-lattice 

interaction  is  still  in  a sense  a small  perturbation,  but  the  last  two  terms 

of  approach  zero.  These  terms  help  bring  about  the  transfer  between 

ep 

external  and  internal  effective  spin  energy  which  is  required  by  the  assump- 


tion of  a canonical  distribution  of  spin  states  with  respect  to 


If  the 


transition  probabilities  due  to  these  terms  are  smaller  than  T-^  , then 


M is  conserved  and  the  treatment  at  the  end  of  the  previous  section 

zp 

applies.  The  fact  that  M is  conserved  in  the  limit  of  small  H.  and  finite 

z p 1 
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H means  that  M = M - M in  this  case;  i.e.,  the  nuclear  magnet- 
er  ZjO  z o ° 

ization  approaches  its  thermal  equilibrium  value,  as  is  expected  on  ele- 
mentary grounds . 

We  may  summarize  these  limitations  with  the  statement  that  the 
theory  is  valid  in  the  range  of  Hj  well  above  the  level  where  the  absorp- 
tion begins  to  saturate. 

Some  of  the  predictions  of  the  theory  are  summarized  in  Fig.  6,  which 

shows  X"  and  8X.'/8H  at  resonance  as  a function  of  H, . At  low  fields  X," 

o l 

is  determined  by  the  line  shape  and  the  Kronig-Kramers  relation  for  zero 
frequency:  oo 

x 

° u v 

so  that 

%M(vo)  = *-  vo^og(vo) 


(44) 

(45) 


The  low  rf  field  dispersion  derivative,  8?0  / 8H  is  also  determined  by 

19  ° 

the  Kronig-Kramers  relations  and  the  line  shape.  We  assume  that  the 

line  is  gaussian.  Using  Eq.(ll)  of  reference  19  (with  the  factory  mentioned 

19  c 

in  the  erratum  ) and  comparing  Eq . (10),  reference  19,  with  Eq . (15), 

reference  25,  we  get  for  a gaussian  line  of  second  moment  <J*H 

8X'/8Ho  = X0Hq/2  <AH2>av  (46) 


This  is  the  low  rf  field  limit  assumed  in  Fig.  6 for  the  case  (curve  a)  of 
pure  dipolar  broadening  (A^=  0).  If  there  exists  an  exchange -type  inter- 
action (A^  ^ 0)  the  resonance  line  will  be  exchange -narrowed  and  the  dis- 
persion derivative  will  be  increased  (curve  b). 


The  region  indicated  by  the  dotted  lines  in  Fig.  6 cannot  be  treated 
theoretically,  but  2C”  and  8%' / 2Hq  in  this  region  are  expected  to  undergo 
a smooth  transition  between  their  low  and  high  rf  field  values,  as  indicated. 


The  solid  lines  in  Fig.  6 for  large  H,  represent  the  predictions  of  (41) 

1 2 
and  (43).  In  the  limit  of  large  the  dispersion  derivative  is  Mq/2H^ 

(assuming  T^  = T^g)  independent  of  the  spin-spin  interaction.  At  inter- 
mediate values  of  H^,  for  pure  dipolar  interaction  (X^  = 0),  (31)  and  (41) 
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predict  that  QX} /BHq  should  approach  a plateau  of  value  3XqHo/4 

(curve  a).  The  presence  of  exchange-type  interaction  will  increase 
2 

6H  and  thus  lower  the  height  of  the  plateau  (curve  b). 

In  the  theory  outlined  above  it  is  implicitly  assumed  that  parts  of  the 
spin-lattice  interaction  can  be  included  in  the  spin  Hamiltonian  as  classical 
perturbations.  These  parts  are  the  Knight  shift  in  metals  and  the  chemical 
shift  in  insulators,  accounted  for  by  replacing  Hq  by  Hq  + AH,  and  the  possible 
mutual  nuclear  spin  interactions  due  to  the  electrons  in  the  solid,  of  which 
the  last  term  of  (12)  is  an  example.  We  know  of  no  rigorous  justification  for 
this  splitting  up  of  the  spin-lattice  interaction  into  stationary  and  relaxation 
parts,  but  it  seems  physically  quite  reasonable. 

33 

In  insulating  crystals,  where  spin  diffusion  usually  plays  an  important 
role  in  the  relaxation  process,  the  theory  above  may  not  be  directly  applicable, 
although  the  qualitative  conclusions  are  apparently  correct.  Spin  diffusion  is 
expected  to  be  affected  by  the  presence  of  the  large  rf  field.  More  specifically, 
(34)  probably  respresents  an  oversimplification  of  the  relaxation  process. 

In  metals,  (34)  may  also  be  an  oversimplification  of  the  actual  relaxation 
process.  Equation  (34)  would  almost  certainly  hold  if  it  were  also  time- 
averaged  in  a suitable  way,  but  it  probably  does  not  represent  the  details 
of  the  nuclear  relaxation  correctly.  In  particular,  the  relaxation  of 
neighboring  nuclear  spins  by  the  conduction  electrons  may  not  be  in- 
dependent, but  may  instead  be  correlated  in  some  way,  owing  to  the  finite 
extent  of  the  electron  wave-functions.  In  this  case  the  reasoning  behind 
(37)  would  probably  still  be  correct,  but  (36)  would  have  to  be  replaced  by 
the  less  specific  equation 

[9/8t)SLR<^er>  = -<^er>/Tlp  . (47) 

where  T^  is  greater  than  g-T  ^ (for  T^e  = T^)  for  small  and  approaches 
T^  (or  in  general  T^jfor  largeH^.  Equations  (30), (32),  (37),  and  (47)  yield 

T , M cos  *9 
M =-lP-°  _ 

°p  Tx(l  + 6H^/Her  ) (48) 


and  at  resonance 
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W = TiP^oHo 
9Hq  2T^if+6i^) 


T,  can  be  experimentally  determined  by  observing  the  phase  of  the 
1 P 

dispersion  signal  relative  to  the  magnetic  field  modulation,  as  discussed  in 
Section  111(e)  below. 

Equation  (47)  may  appear  at  first  sight  to  be  inconsistent  with  (36),  when 

the  spin  relaxation  is  incoherent.  Actually,  this  is  not  the  case,  since^^gg)^ 

and  M H are  related  to » by  (30)  and  (32).  If  we  solve  these  four 
z p er  Ler  ' 

equations  for  T j we  get  the  same  expression  obtained  in  Section  111(e)  in 
connection  with  the  theory  of  fast  modulation  with  incoherent  relaxation 
(Eq.  63). 

32 

Bloembergen  has  suggested  a method  of  calculating  T in  those 

metals  for  which  the  nuclear  spin-^pin  interaction  A..  I.*  I.  is  the  pre- 

-ir  26  27  JK~3 

dominant  term  in^^g.  This  is  the  case  ’ in  most  of  the  heavier 

metals,  probably  including  copper  and  aluminum.  Equation  (47)  can  be 

written  (assuming  Tj  = ^2e^ 

IjT  W<?4r»  = Mz„Her/T  1 ' <>fSS>  ^SS  ' <50> 

The  term  MzpHer /T  j in  (50)  is  deduced  in  the  same  way  as  in  ( 36) . . Tgg  * 
is  a suitable  average  of  the  probability  that,  of  two  neighboring  nuclei, 
either  one  will  be  flipped  by  an  electron  without  the  other.  If  the  electron 
wavelength  is  very  short  it  is  expected  that  the  relaxation  will  be  incoherent 
and  that  Tsg  = iTj  as  in  ( 36) . 

T can  be  estimated  by  considering  nearest  neighbors  separated  by  a 
s s 

distance  The  important  nuclear-electronic  interaction  is  A[ (J.  6(jj-^r. ) 

+ I,  6(r-r,  )1  *S  , where  S is  the  electronic  spin  operator  and  A is  an  inter- 

/w  lr  '*-v  <vv  j /W\  ' KM  1 1 

action  constant.  For  a transition  involving  k and  k'  = k + Ak  as  initial 
and  final  electronic  propagation  vectors  this  interaction  becomes 


APp(e 


J I.  + e 
/v'0 


i Ak*  r 
~ «k 


= 2 e 


apf  +ik)[cos  + 1 sin  ^'Xjk* + nj- 


(51) 
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+ (I.  - I,  )[  cos  Ak*  r +i  sin  { Ak-  r.i  ) - 1] 

VAj  ~vk,L  *v~vjk  AVASjk'  ■> 


(51) 


t.  + I,  commutes  with  A.,  I. -I,  and  thus  does  not  contribute  much  to  the 
Wj  ^ 'UL 

decay  of  <</Tsg»  unless  the  term  Bilrl.^iw  is  relatively  large.  Assuming 


jk  jz  kz 


B^k  small,  most  of  the  relaxation  results  from  the  term  involving 

~j  " Jk’  transition  probabilities  proportional  to  2-2  cos  We 

assume  that  the  energy  contours  in  k-space  are  spherical,  in  which  case 

k — k1  ^k  , the  Fermi  level  value  of  k,  and  that  all  directions  in  k-space 
m 

are  equally  probable  for  the  initial  and  final  states  . is  then  proportional 

to  the  average  of  2-2  cos  (k-k'fr^  over  all  directions  of  k and  k' , keeping 


k = k'  = k , and  is  normalized  by  the  requirement  T, 


SS' 


-*T! 


as  k 


m 


In  this  way  we  get 


r"1 

ss 


— 2T,1  [ 1-  ( 1 - cos  2 k r. . ) / 2k^  r ?. ] 
1 1 ' m ij  m ijJ 


CD  . 


(52) 


If  it  is  assumed  that  there  is  one  electron  per  atom  in  the  Fermi  gas, 

then  k r.,  = 3.  36  for  a face-centered-cubic  crystal  and  (52)  indicates  that 
m jk  ’ 

Tgs~ J^i*  yielding  almost  the  same  predictions  as  (41)  and  (43).  Only  if 

k rM  <tr/2,  or  if  the  conduction  band  is  almost  filled  and  can  be  described 
m jk 

as  an  almost  empty  band  of  positive  holes  with  k r..  < tt/2,  will  there  be 

m j k 

an  appreciable  effect  due  to  coherent  relaxation  of  neighboring  nuclei. 


D.  Two  Magnetic  Ihgredients 

There  are  two  limiting  cases  to  be  considered  when  there  are  two  or 
more  types  of  spin  in  the  lattice. 

The  first  of  these  is  typified  by  F centers  whose  electron  spin  resonance 
is  broadened  by  hyperfine  interaction  with  the  surrounding  nuclear  spins. 
Although  the  different  F centers  in  an  alkali  halide  are  magnetically  coupled 
in  theory  via  nuclear  spin  diffusion  and  direct  interaction  with  each  other, 
in  fact  such  coupling  is  negligible  over  a length  of  time  comparable  to  T j 
(electronic)  because  the  spin-spin  transitions  involved  are  highly  forbidden 
energetically  by  the  relatively  large  local  field  differences  experienced  by 
the  different  spins.  Then  the  F-centers  must  be  regarded  as  decoupled 
and  the  assumption  of  a canonical  distribution  of  states  with  respect  to 

n 

the  transformed  Hamiltonian  abandoned.  Portis  has  developed  an  exact 
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theory  for  this  case. 

The  other  limiting  case  is  typified  by  two  different  nuclear  species  in 

the  same  crystal.  In  this  case  the  local  field  differences  between  nuclei  are 

relatively  small  and  spin  diffusion  can  take  place.  Following  Van  Vleck  we 

ufe  primes  to  denote  the  spins  whose  resonance  is  not  being  observed;  the 

rf  frequency  to  is  supposed  to  be  near  the  resonant  frequency  of  the  unprimed 

25 

spins.  The  spin  Hamiltonian  is  given  by  Van  Vleck.  We  transform  to  a 
rotating  system  as  before,  using  the  transformation  t in  which  the  sum 
over  the  1.^  is  carried  out  over  all  the  nuclei,  primed  and  unprimed.  The 
transformed  Hamiltonian  is 

Xev  - Wfler'  + g<AjkV4  + VjAz>  +SCjKyk,a 


(53) 


Here  g'  and  y'  are  respectively  the  g factor  and  gyromagnetic  ratio  of 

the  primed  ingredient,  A.,.,  and  B.,,  , are  given  by  (22)  and  (23)  and 

J K J K 


Cjk-  = xjk*  + (1  * 3^jk‘  > ss’P 


2 -3 
r Jk’  J 


(54) 


In  (53)  we  have  neglected  time-de pendent  terms  as  before.  We  can  also 

neglect  the  term  g 1 PH ^ which  is  nonsecular  if  | w-y'Hj  is  much 

greater  than  the  resonance  line  width -of  the  unprimed  nuclei.  The  term 

involving  2Zl.>  then  commutes  with  the  rest  of  the  Hamiltonian  and 
J J* 

can  be  ignored,  since  the  unprimed  spins  are  not  affected  directly  by 
the  rf  field. 

The  remaining  terms  of  the  Hamiltonian  are  secular  perturbations  and 
do  not  commute  with  each  other.  This  means  that  they  are  coupled  to- 
gether, with  the  rather  surprising  result  that  effective  energy  can  be 
transferred  from  the  unprimed  system  to  the  primed  system.  When  rf 
energy  is  absorbed  by  the  unprimed  spins,  part  of  this  energy  will  be 

transferred  to  the  primed  system  via  the  interaction  /,  . C ,1.  I , 

j ,k'  JK-  JzK 

This  energy  will  be  entirely  in  the  form  of  spin-spin  (internal)  energy; 
the  external  energy  of  the  primed  spins  (expectation  value  of  the  last 
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term  of  yC  ) will  be  unchanged.  The  average  populations  of  the  primed 
uer 

spin  levels  will  remain  at  their  equilibrium  values,  but  the  primed  spins 
will  become  more  ordered  in  orientation  with  respect  to  their  local  fields. 


The  statements  in  the  previous  paragraph  may  appear  to  contradict 
the  usual  assumption  that  two  different  nuclear  magnetic  ingredients  in  a 
solid  interact  entirely  independently  with  the  rf  field,  from  which  it  would 
follow  that  the  primed  system  is  unaffected  by  rf  power  at  the  resonant 
frequency  of  the  unprimed  system.  Actually  this  assumption  is  not  quite 
true,  since  the  spin-spin  energy  of  the  primed  system  evidently  increases 
in  this  case.  For  most  purposes  this  energy  is  negligibly  small  compared 
to  the  external  energy  (in  the  fixed  system)  of  the  primed  system,  which  is 
unaffected  by  the  rf  field.  Only  in  the  present  case,  where  effective  energy 
is  the  important  quantity,  will  the  spin-spin  energy  be  important,  and  then 
only  in  determining  the  behavior  of  the  unprimed  nuclei. 

As  before  we  assume  that  at  high  rf  levels  the  spin  system  is  in  a 
canonical  distribution  of  states  with  respect  to  the  retained  parts  (all  but 
the  last  term)  of  the  effectively  time-independent  transformed  Hamiltonian 
^=('er-  As  a consequence  of  this  assumption,  the  part  of  the  nuclear 
magnetization  due  to  the  unprimed  spins  is  in  the  effective  field  direction 
with  magnitude  MZp  while  that  due  to  the  primed  spins  remains  at  its 
thermal  equilibrium  value  in  the  z-direction. 


The  effect  of  the  unlike  neighbors  can  be  expressed  in  terms  of  the 
2 2 

ratio  6H  /Ho.  of  internal  spin-spin  energy  to  external  (unprimed)  spin 


er 


energy.  In  Appendix  C it  is  shown  that 


6H2  = J<AH2>u  + I(I+l)(fNh2y2)_1(^A2  +X1  aL)  + <AH2> 

j*k  J jfck'  J * 


+ ^-(f/f')<4H2>_ 


(55) 


where <^AH  ^ is  the  contribution  of  the  unprimed  ingredient  to  the  second 
25  u 2 

moment  of  its  resonance,  ^^1  \ is  the  contribution  of  the  primed  mag* 

'pn  2v 

netic  ingredient  to  the  second  moment  of  the  unprimed  resonance,  } 

is  the  contribution  of  the  primed  ingredient  to  the  second  moment  of  the 

primed  resonance,  N is  the  total  number  of  both  ingredients,  and  fN  and 
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f'N  are  the  jiumbers  of  unprimed  and  primed  nuclei  respectively. 

If  the  relaxation  times  of  the  two  types  of  nuclei  are  equal  the  quasi - 
static  equilibrium  value  of  the  unprimed  magnetization  is  given 

by  (38)  using  the  same  reasoning  as  before,  with  taken  to  be  the  unprimed 
contribution  to  the  static  susceptibility.  If  the  relaxation  times  of  the  un- 
primed and  primed  nuclei  are  and  T|,  then  it  can  be  shown  (Appendix  C) 
that 

M = M cos®  -fl+H-2  [26H2  +26H2  T./T'  +<2H2>  (l+T./T')]!" 
opo  I er  L u pl.x  pu  1 1 ( 

1 ' J (56) 

2 2 2 
where  6 is  equal  to  the  first  two  terms  of  6 H , (Eq.  55)  equals 

the  third  and  fifth  terms  of  6H2,  and  is  the  unprimed  contribution  to 

the  static  equilibrium  magnetization. 

As  before,  the  theory  is  valid  only  for  rf  field  levels  above  saturation, 
and  if  the  spin-lattice  interaction  results  in  correlated  scattering  for 
neighboring  spins  the  remarks  at  the  end  of  Section  III(c)  apply. 

The  predicted  behavior  of  X”  and  93C/9Hq  at  resonance  for  a crystal 
with  two  magnetic  ingredients  is  shown  in  Fig.  6.  The  absorption  is  as 
before,  except  that  g(vQ)  may  be  different  and  the  parameters  of  the  satu- 
ration curve  refer  only  to  the  unprimed  nuclei.  The  predicted  dispersion 
shown  is  that  for  a crystal  in  which  the  observed  (unprimed)  nuclei  are 
in  the  minority,  the  relaxation  times  T,  and  TJ  are  comparable,  and  the 
absorption  is  gaussian  below  saturation.  The  plateau  in  8%'/SHo  above 
saturation  is  reduced  relative  to  that  for  the  pure  case  (curve  a)  because 
of  the  additional  spin-spin  interaction  of  the  unlike  nuclei.  The  line  shape 

above  saturation  is  the  same  as  for  one  magnetic  ingredient,  with  a suitable 
2 

increase  in  6H  . 

E . Fast  and  Intermediate  Modulation 

When  the  period  of  the  magnetic  field  modulation  becomes  comparable 
to  the  spin-lattice  relaxation  time  the  observed  absorption  and  dispersion 
signals  above  saturation  will  differ  from  those  obtained  if  T ^ is  very  short. 
The  same  remark  applies  if  the  rf  amplitude  or  frequency  is  modulated.  In 
solids  it  is  difficult  to  avoid  modulation  effects  because,  except  for  the 


r 

- — 
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heavier  metals,  the  nuclear  relaxation  time  at  300  K.  is  greater  than  10 
milliseconds,  requiring  modulation  frequencies  co^  of  10  cycles  or  less. 

Such  low  modulation  frequencies  are  seldom  used  in  practice,  because  of 
noise  and  stability  considerations.  This  is  evidently  why  the  non-saturation 
of  the  dispersion  reported  here  has  not  been  previously  observed. 

Another  reason  for  considering  fast  modulation  is  that  by  so  doing  it 

is  possible  to  obtain  the  spin  lattice  relaxation  time.  This  possibility 

16 

has  been  exploited  by  Halbach,  ’who  has  developed  the  theory  of  fast 
magnetic  field  modulation  for  systems  obeying  the  Bloch  equations. 

5 6 

Theories  ’ of  fast  modulation  based  on  the  assumption  of  a spin 

temperature  lead  to  incorrect  results,  even  for  systems  obeying  the  Bloch 

equations.  In  these  theories  it  is  assumed  that,  for  the  spin 

temperature  is  that  corresponding  to  the  average  values  of  H , H,,  andu  . 

- 1 o i 

Actually,  for  T j yH  ^ the  nuclear  magnetization  will  remain  in  the 

direction  of  the  effective  field,  and  thus  its  z-component  and  temperature 

will  vary  sinusoidally  with  the  z-component  of  the  effective  field. 

In  this  section  we  treat  only  modulation  of  the  magnetic  field  Hq. 
Frequency  modulation  of  the  rf  field  is  equivalent,  and  amplitude 
modulation  of  Hj  can  be  treated  similarly.  We  assume  that  to^^yH^  and 
that  u is  much  smaller  than  the  unsaturated  line  width,  so  that  modulation 


m 


34 


effects  of  the  type  considered  by  Kar plus’'*  are  negligible. 

The  applied  magnetic  field  is 

H = (H  + H cos  oo  t)z  + (2H,  cos  cot)x  • 
o m m~-  1 ~r 


(57) 


In  addition  to  the  explicit  time  dependence  of  H,  in  (57), 'to  also  increases  or 
decreases  very  slowly  in  time  as  the  spectrometer  sweeps  through  the 
resonance  line. 

In  general  the  ^component  of  magnetization  is  (to  fir  st  order) 

M = (M1  + M'  , costo  t + M'  , sinu  t)  coswt 
x ' x xl  m xZ  m 

(5«) 

+ (MM  + M" . cos  to  t + M",  sin  w t)  sinut  > 

' x xl  m x2  m 

where  M^.,  M^,  etc.  are  constant  in  time  (or  vary  only  very  slowly  as  to  is 
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varied) . 

The  apparent  dispersion  derivatives  are  defined  as 


a*'/3Ho  = . 

(59a) 

sxymo  = Mi2/2HlHm  , 

(59b) 

(59c) 

= M^2/2HlHm  . 

(59d) 

8%,' /8H  and  8%''/3H  are  the  apparent  dispersion  and  absorption  derivatives 
1 o lo 

observed  when  the  lock-in  detector  is  adjusted  in  phase  with  the  modulation, 
and  8X^/8Hq  and  8%'_j/8Ho  are  the  apparent  dispersion  and  absorption 
derivatives  observed  when  the  lock-in  detector  is  adjusted  90°  out  of  phase 
(in  quadrature)  with  the  modulation. 

The  Hamiltonian  is  transformed  to  a rotating  system  as  before,  but  ® 
and  Her  depend  on  time  sinusoidally  with  frequency  wm-  This  frequency  is 
so  low  that  it  induces  no  appreciable  transitions  between  the  eigenstates  of 
, and  it  can  be  regarded  as  a reversible  adiabatic  (slow)  perturbation 
of  the  spin  system  (not  the  whole  system,  as  1).  It  is  then 

reasonable  to  assume  as  before  that  the  spin  system  remains  in  its 
highest  entropy  macrostate;  i.e.,in  a canonical  distribution  of  states. 

In  Appendix  B it  is  shown  that  this  assumption  leads  to  the  equation 
(for  one  magnetic  ingredient) 


4« 


ar 


■«H?>  - 


M H /T,  - Z«J4CJ>  /T,  - M cos0H  /T, 
zperl  SS  1 o er  1 


-M  -H 
/w\ x ~*>er 


(60) 


Using  (32)  gives 

= (+M  H - M cos  ® H - 2<<>4C.»)/T1 +M  -H  (61) 

dtSS  zper  o er  ^ SS  1 /~r  ««er 

35 

Equation  (61)  also  follows  directly  from  the  first  law  of  thermodynamics 

for  the  transformed  Hamiltonian  dU1  = dQ  + H-  dM,  where  U'  is  the 

er  er  w er 

effective  internal  energy  of  the  spin  sy stem^J^g^) , and  Qer  is  the  effective 
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heat  transferred  from  the  lattice  to  the  spin  system  via  the  spin-lattice 
interaction.  Application  of  thermodynamics  to  is  valid  because  the 

changes  in  are  reversible  (slow)  and  the  coupling  to  the  lattice  is 

assumed  weak. 


The  assumption  of  a canonical  distribution  of  states  implies  as  before 
that  the  magnetization  is  in  the  direction  of  i}er  with  magnitude  n*Zp>  and 
that<^Jk£ec»  = -M  (6H2/H  ).  Substituting  into  (61)  we  get 

bb  IS  CJ*  ^ 


where 


M = -(M  -M  )/T,  +M  H H"1  (1 +H2  / 5H2) 
zp  ' zp  op  lp  zp  er  er ' er  v 


Tlp  = Tx(l  + 6H2/H2r)(l  + 26H2/H2r) 


-1 


(62) 


(63) 


and  is  given  by  (38). 

Similar  reasoning  shows  that  when  two  magnetic  ingredients  are 
present  (62)  still  holds,  and  Mq  and  are  given  by  (48)  and  (56).  If 
the  spin-lattice  relaxation  is  correlated  between  neighbors  the  considerations 
at  the  end  of  Section  III(c)  apply,  and  is  the  same  quantity  introduced 

there.  M ^is  again  given  by  (48)  and  T^p  can  be  empirically  determined 
as  described  below. 

If  T^p  is  large  compared  to  the  time  taken  to  sweep  through  the  entire 
resonance,  we  can  neglect  the  term  (M^p  - M )/T^pin(62).  Assuming 
that  initially  w is  well  off  resonance  (|Hq  -w  /yp»6H,  Hx)  and  + Mz  = Mz=Mo 
we  get 


+ M 


M 


zp 


TT 


(1+  6H7Her)I 


(64) 


and 


9X’  / 9H  = + 9/9H 
i o - o 


MQsin®/  2H1 


(1+ 


sriTH2  )i 


er 


(65) 


9X! L/9H  = dX'J  / 9 H = 9%"/9H  = 0. 

2 o 1 o Z o 


(66) 


The  plus  or  minus  sign  in  (64)  and  (65)  depends  on  whether  u>  is  initially 
less  or  greater  than  the  resonance  frequency.  This  is  the  case  Bloch^ 
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calls  adiabatic  fast  passage.  The  present  analysis  shows  that  to  observe 

a fast  passage  siganal  it  is  not  necessary  to  sweep  through  the  resonance  in 

4 

a time  short  compared  to  T^>  as  previously  supposed,  but  only  in  time 
short  compared  to  T^.  In  the  case  of  solids,  where  T^»T^,  this  makes 
observation  of  fast  passage  signals  much  easier.  It  is,  of  course,  still 
necessary  to  use  rf  levels  well  above  the  saturation  level. 

If  T . = T , is  short  compared  to  the  time  taken  to  pass  through  resonance. 

Ip  l 

we  can  regard  co  as  fixed,  and  the  first-order  solution  of  (62)  is 


M = Kt  + M , cosu  t + M,  sin« t,  (67) 

zpop  1 m Z m 

where  ^Qp  i8  the  quasi -steady-state  value  of  = M corre8Pon(^n6 

to  Hq  = Hq.  and  are  given  by 


, , - 1 M H cos  © 

M,=  — (1+liti  T 7 ) ( °p  5 -r-~, 

1 m lp  H (1+H2  / 6H2)  o>  ZT2,  8H 

e r e r ml p 


M p = L-t-M.  + H 

2 “mTlp  1 


3M 


m 


■SH- 


OP 


) • 


H 3M  _ 
m op. 


(68) 

(69) 


Here  3M  /3H  is  evaluated  for  H = H . (68)  and  (69)  are  obtained  by 

op  o o o 

substituting  (67)  into  (62)  and  setting  = (Hmcosumt)3M^  /3Ho+Mo^. 

The  apparent  dispersion  derivatives  are  (from  Eqs  . (33),  (56),  (59), 
and  (67)  ) 

dX / 3H  = M',  sin  ®/2H,H  -M  sin2  © cos  @/2H?,  (70) 

lol  l m op  i 

QX'/dH  = ML  sin  ®/2H.H  . (71) 

Z o 2 1 m 


At  resonance. 


- 1 


dX\IBHo  = (1  +“mT  lp>  «'/8H0, 

8;l;2'8Ho=“mTlp  M<V8Ho  . 


(72) 

(73) 


Here  9X*  / OH  is  the  true  dispersion  derivative  at  resonance  in  the  limit 
o 

of  slow  modulation.  Equation  (7  3)  shows  that  the  relaxation  time  T, 


can 
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be  experimentally  determined  from  the  ratio  of  the  in-phase  to  the 

quadrature  dispersion  signals,  and  the  true  dispersion  derivative  can 

be  calculated  from  (72).  This  method  of  measuring  the  relaxation  time 

18 

was  previously  suggested  and  applied  to  liquids  by  Halbach.  Similar 

36 

effects  have  evidently  also  been  observed  by  Portis  and  Shaltiel.  It 
will  be  noted  that  the  relaxation  times  can  also  be  measured  by  modulating 
the  rf  power  level,  as  is  sometimes  more  convenient  for  microwave  studies 
of  paramagnetic  relaxation. 


The  apparent  absorption  derivatives  are  obtained  from  conservation 
of  energy  considerations  in  the  fixed  coordinate  system.  The  rate  of 
transfer  of  energy  from  the  spins  to  the  lattice  is  to  a good  approximation 
H (M  -M  )/T,  . Energy  is  absorbed  by  the  spins  from  the  external  field 

° ° ■ . 1 //-J  y 

at  a rate  H*M.  The  internal  energy’Wtgj,^  also  varies  with  the  modulation 
but  this  variation  is  of  the  order  of  6H2/HQHer  times  the  external  energy 
change  II* hi  and  can  be  neglected.  Conservation  of  energy  then  requires 
that 


+H-M  = H (M  -M  )/T,  . 
o'  o Z 1 


A/V\ 


(74) 


Averaging  over  one  rf  period  we  get 

H (M  -M  ) /T  , = - uH.(M"  + M"  , cos  w t + M"  sin  w t)  + H hi, 
o'  o z 1 1'  x xl  m x2  m ’ o 2 

(75) 

Equations  (33),  (58),  (59)  and  (75)  can  be  solved  for  the  apparent  absorption 
derivatives.  The  result  is  complicated  and  will  be  omitted.  The  ab- 
sorption signal  depends  on  both  the  modulation  frequency  and  phase.  Since 
it  is  not  customary  to  adjust  a lock-in  detector  with  extreme  care,  the 
phase  of  the  observed  signal  in  previous  measurements  of  T ^ by  saturation 
must  be  regarded  as  uncertain,  and  the  reported  values  of  T j correspondingly 
uncertain.  The  resulting  errors  in  T j are  not  likely  to  be  greater  than  a 
factor  of  two  or  three  because  the  onset  of  saturation  will  still  occur  at 
level  H1^2[y2T1g(vQ)]‘  2 corresponding  to  the  point  where  the  rf  transition 
probability  is  comparable  to  the  spin-lattice  transition  probability.  Fast 
modulation  will  only  change  the  details  of  the  absorption  saturation  and 
apparent  asymptotic  absorption  at  resonance  for  large  H^. 
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The  predicted  dependence  of  Tj^  on  rf  field  is  shown  in  Fig.  6, 
assuming  uncorrelated  spin-lattice  relaxation  between  neighbors.  In  the 
limit  of  large  Hj,  T^  equals  T ^ (or  more  generally  here  we  assume 

as  usual  At  the  value  of  corresponding  to  the  knee  of  the 

dispersion  curve  undergoes  a transition  to  T^/2  for  a single  magnetic 

ingredient,  or  to  some  other  value  for  two  magnetic  ingredients  as  indicated 
by  (48)  and  (56).  Below  saturation  the  observed  dispersion  signal  is  ex- 
pected to  be  in  phase  with  the  modulation  and  to  correspond  to  the  true 
dispersion  derivative.  The  observed  T ^ as  defined  by  (73)  is  then  ex- 
pected to  decrease  in  some  unpredicted  way  corresponding  to  the  dotted 
lines  in  Fig  . 6 . 


IV.  Comparison  with  Experiment 


We  have  not  made  a detailed  analysis  of  the  observed  line  shapes,  but 

it  appears  that  they  are  in  agreement  with  the  predictions  of  the  theory. 

In  the  metals  the  absorption  and  dispersion  derivatives  above  saturation 

19 

are  very  nearly  Lorentzian  and  have  the  expected  width.  The  apparent 
dispersion  derivatives  in  NaCl  at  large  rf  fields  appear  to  be  in  agreement 
with  the  theory,  assuming  umTj»l  and,  for  the  impure  NaCl,  assuming  T j 
short  compared  to  the  few  minutes  taken  to  sweep  through  the  resonance, 
or,  for  the  Harshaw  NaCl,  assuming  T ^ comparable  to  the  time  taken  to 
sweep  through  the  resonance. 


The  observation  of  the  adiabatic  fast  passage  signal  in  the  Harshaw 
NaCl  is  evidence  that  we  were  justified  in  neglecting  the  time-dependent 
terms  in ^eT  and  (Eqs . (20)  and  (21)  ).  If  these  terms  could  induce 

transitions  among  the  different  eigenstates  of  the  time-independent  part 
of  the  transformed  Hamiltonian  then  the  magnetization  in  the  z ^direction, 

Mz  p,  would  be  destroyed  when  w passed  through  resonance.  The  observations 
onthe  Harshaw  NaCl  indicate  that  the  relaxation  time  of  Mzp  produced  by 
these  terms  is  greater  than  about  one  minute,  so  that  they  can  be  neglected 
for  most  purposes. 


We  now  reconsider  the  absorption  data  of  Fig.  1.  In  the  limit  of  large 

2 

H^,  (43)  indicates  that  X"  should  approach  %oHq/2T  ^yH^  . In  the  limit  of 
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small  Hj  the  Kronig-Kramers  relations  indicate  that%"  should  approach 

1 /2tr JCqv Qg( v),  where  g(v)  is  the  normalized  unsaturated  line  shape.  An 

attempt  has  been  made  to  draw  these  asymptotes  for  the  data  of  Fig.  1. 

2 2 

The  asymptotes  should  cross  at  the  value  of  given  by  l/2-y  HjT^g(vQ)  =1. 
In  this  way  we  get  5.5  milliseconds  for  T ^ aluminum  and  3.55  for  copper. 
Since  it  is  uncertain  whether  the  asymptotes  drawn,  in  Fig.  1 are  the  true 
ones  these  values  of  must  be  regarded  as  upper  limits.  In  addition, 
inaccuracies  in  the  absorption  data  may  introduce  as  much  as  20  per 
cent  error  in  T ^ . These  values  of  T ^ are  not  necessarily  more  reliable 
than  those  obtained  in  Section  II,  and  do  not  appreciably  alter  the  con- 
clusions reached  there  concerning  the  electronic  structure.  They  lead 

to  values  of  v P„  rrJ'/m  of  220  for  aluminum  and  240  for  copper, 
o F 

The  dispersion  data  of  Figs.  2 and  3,  together  with  (72)  and  (7  3), 
yield  values  of  the  true  dispersion  derivative  contribution  per  nucleus 
shown  in  Fig.  7.  Also  shown  in  Fig.  7 is  the  relative  dispersion  de- 
rivative per  proton  for  protons  in  water  doped  with  paramagnetic  im- 
purity. This  is  the  same  H,0  line  shown  in  Figs.  2 and  3.  To  get  the 

- 1 ^ 

absolute  magnitude  of  N 8%V9Hq  for  comparison  with  experiment  we  can 
use  either  the  I^O  data  or  the  aluminum  or  copper  dispersion  derivative 
below  saturation  as  a calibration. 


In  order  to  use  the  proton  dispersion  as  a calibration  it  is  necessary 

22 

to  assume  a value  of  ^H^IH  (5)»  Wangsness  and  Bloch  predict 

that  this  ratio  is  unity  for  protons  in  water  relaxed  by  paramagnetic  im- 

32  37 

purities,  as  considered  here . The  same  result  follows  ’ from  a 

23  37 

somewhat  more  general  theory  of  Kubo  and  Tomita,  Solomon  has 
recently  measured  and  T using  spin  echo  techniques  for  protons 

in  water  containing  Fe  ions.  For  Fe  concentrations  up  to  that  required  to 
reduce  Tj  to  one  millisecond  he  finds  that  = + *03,  in 

agreement  with  theory.  Unfortunately  the  dispersion  derivative  observed 
in  water  in  Fig.  7 is  inconsistent  with  that  observed  below  saturation  in 
copper  and  aluminum  (which  is  rigorously  determined  by  the  Kronig-Kramers 
relations  and  the  unsaturated  absorption  line  shape)  unless  we  assume 
1H^T  2H~^ ' Thus  it  appears  that  either  the  errors  in  the  data  are 
greater  than  the  conservative  estimate  of  + 20  per  cent,  or  that  there  is 


r 
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something  wrong  with  the  theory  of  nuclear  resonance  in  liquids.  The 

latter  possibility  should  not  be  taken  too  seriously  in  view  of  the  rather 

preliminary  nature  of  the  data,  but  it  may  be  that  the  high  concentration 

of  paramagnetic  impurites  (—1  per  cent)  could  cause  a decrease  in  the 

intensity  of  the  dispersion  signal  because  the  r.m.s.  value  of  the  rapidly 

fluctuating  local  fields  seen  by  the  nuclei  is  larger  than  the  applied  field 

H . 
o 

In  view  of  this  difficulty  we  adopt  the  less  accurate  (but  theoretically 
more  rigorous)  alternative  of  using  the  dispersion  signal  in  aluminum  and 
copper  below  saturation  as  a calibration.  We  approximate  the  line  shape 
with  an  inverted  truncated  parabola: 


g(v)  = 3a  ^ (aZ  - vZ)/2  for  vZ<aZ 


(76) 


g(v)  = 0 


, 2 2 
for  v >a 


This  is  an  excellent  approximation  for  aluminum  and  a fair  one  for  copper, 
which  has  a more  nearly  gaussian  resonance  line  than  aluminum.  Inserting 
(76)  into  the  Kronig  -Kramers  relations  and  evaluating  the  second  moment 
leads  to 


or 


N 


.-1 


W-  = -«0Ho'<AHJ>av 


dX' 


= . 6ir  I(I+l)v  y/3kT  <AHZ> 


AV 


(77) 


(78) 


where  ^AH  is  the  second  moment  of  the  line.  If  the  line  were 

assumed  to  be  gaussian  the  factor  in  (7  8)  would  be  1.0  instead  of  .6.  In 
copper  (78)  predicts  a dispersion  derivative  which  is  expected  to  be  slightly 
too  small,  because  of  the  more  nearly'gaussian  shape  of  the  copper  reso- 
nance . 

2 2 

In  aluminum  the  predicted  second  moment  <^AH  ^ is  7.5  gauss,  . 

— ' A V 


26 


neglecting  possible  anisotropic  electronic  coupling  between  nuclei. 

1 2 
Gutowski  and  McGarvey  report  an  experimental  value  of  10.5  gauss 

39  2 2 

and  Rowland  a value  of  8.7  gauss.  We  assume  a value  of  9 gauss. 
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In  copper  Gutowski  and  McGarvey  report  that  experimentally 

is  6.3  gauss  / in  reasonable  agreement  with  their  theoretical  prediction 

2 1 
of  5 . 6 gauss.  We  use  their  experimental  value.  Corcrections  to 

^^y  which  depend  on  Tj  are  negligible  compared  to  the  other 
approximations  and  uncertainties  in  the  theory. 

We  use  the  aluminum  dispersion  derivative  as  a primary  calibration, 
and  assume  that  its  experimental  value  (upper  horizontal  arrow  in  Fig.  7 
is  given  correctly  by  (78).  The  copper  dispersion  derivative  below 
saturation  is  then  somewhat  greater  than  the  prediction  of  (7  8)  (lower 
horizontal  arrow),  as  expected.  The  discrepancies  between  the  ex- 
perimental points  for  aluminum  at  small  are  almost  certainly  due 
to  experimental  error  . 


The  dotted  lines  in  Fig.  7 are  the  predicted  high  rf  field  asymptotes 
for  the  dispersion  derivatives,  again  using  the  low  rf  field  aluminum  dis- 
persion as  a calibration.  The  predicted  asymptote  is  N ^SX'/SH  > 

2 ^ 
XHq/2NH1,  as  obtained  from  either  (10)  or  (41),  assuming  Tj  = The 

solid  curves  in  Fig.  7 for  large  H,  correspond  to  the  prediction  of  (41), 

2 2 1 2 
taking  6H  =3.4  gauss  for  aluminum  and  5.0  gauss  for  copper.  In 

aluminum  the  plateau  predicted  for  8X'/9Hq  at  intermediate  fields  is  not 

resolved,  owing  to  the  relatively  small  ratio  of  T ^ to  (y  6H)  1 compared  to 

that  assumed  for  the  predictions  of  Fig.  6.  In  copper  the  plateau  evidently 

nearly  coincides  with  the  value  of  8X1  /8H  below  saturation. 

o 

In  Fig.  8 are  plotted  experimental  values  of  T ^ obtained  from  (7  3) 
and  the  data  of  Figs.  2 and  3.  For  large  the  observed  T ^ values 
agree  fortuitously  well  with  the  values  of  T ^ obtained  in  the  beginning  of 
this  section,  and  disagree  with  the  values  of  T ^ in  Table  1. 


The  behavior  of  T is  in  conflict  with  (63),  which  is  based  on  the 
IP 

assumption  of  incoherent  relaxation  of  neighboring  nuclei  implicit  in  (35). 

If  this  assumption  were  correct  T,.  should  reach  a plateau  of  about  4t. 

1 p 2 1 p 

at  the  value  of  H^(  — 1 gauss)  at  which  the  intermediate  plateau  of  the 
dispersion  derivative  occurs.  Such  a plateau  apparently  exists  but  the 
value  of  T^  is  too  large.  Such  a large  discrepancy  cannot  be  explained 
by  coherent  nuclear  relaxation  unless  the  electronic  band  structure  is 
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appreciably  different  from  the  free  electron  approximation  with  one 
electron  per  atom.  We  do  not  know  the  reason  for  this  behavior. 


Since  the  discrepancy  in  Tjp  is  evidently  real,  at  least  in  the  case 
of  copper,  we  are  forced  to  use  the  more  empirical  approach  described 
at  the  end  of  Section  HI.  For  1 gauss  T ^ can  be  approximated  by 


T2e(l  + 6H2/H2)  (1  + a6H2/H2  ) 


(79) 


where  T,  is  5.5  milliseconds  for  aluminum  and  3.55  milliseconds  for 
copper,  and  a is  1.7  for  aluminum  and  1. 3 for  copper.  T2g  is  used  in 
(79)  because  T.  is  expected  to  approach  T,  , the  transverse  spin-lattice 

1 p ZG 

relaxation  time,  in  the  limit  of  large  Hj  and  at  resonance. 

If  we  assume  that  T^  = T2g  as  would  be  indicated  by  the  values  of 

T.  obtained  at  the  beginning  of  this  section,  (41)  predicts  the  same 

1 2 
depencence  of  dX.'  /9H  as  before  (solid  lines  in  Fig.  7)  if  we  take  6H  =4 

2 ° 2 
gauss  for  aluminum  and  7.7  gauss  for  copper. 


A third  interpretation  of  the  data  obtains  if  we  assume  that  the  values 

of  T.  in  Table  1 are  correct.  This  assumption  is  in  conflict  with  the 
1 12  22 

prediction  ’ Tj  = T2g  for  a nucleus  in  a Fermi  gas  of  electrons.  The 

predicted  asymptote  of  8^' /9Hq  is  increased  by  T^g/T  j and  (41)  is  in 

reasonable  agreement  with  the  data  if  we  take  6h^  = 5.9  gauss2  for 

2 

aluminum  and  10  gauss  for  copper. 


We  now  estimate  the  theoretical  value  of  6H  . If  we  assume  A.,  =0, 

2 2 JK 

6H  is  predicted  by  (31)  and  (55)  to  be  2.5  to  3 gauss  for  aluminum  and 
2 

2.85  to  3.1  gauss  for  copper.  The  lower  limits  of  these  estimates 

correspond  to  the  classical  dipolar  interaction  between  nuclei  and  the 

1 38 

upper  limits  correspond  to  the  fact  that  the  observed  ’ second  moments 

of  the  resonance  line  are  ten  to  twenty  per  cent  larger  than  their  theoretical 

26) 

values,  possibly  because  of  additional  psuedo-dipolar  coupling  between 
nuclei  via  the  conduction  electrons  . 


It  is  clear  that  the  observed  values  of  6H  cannot  be  explained  by 
classical  dipolar  coupling  alone.  Ruderman  and  Kittel  have  calculated 
the  magnitude  of  the  exchange -type  coupling  between  nuclei  due 
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to  the  conduction  electrons;  the  result  of  their  calculation  can  be  written 
(using  Eq.  (2)  ) 


rsj 

A.,  = 
Jk 


v^1i[2k  r..  cos  (2k  r..)-sin(2k  r.,)1 

o 1 m jk  ' m jk'  ' m jk,J 

"Z 


Sir  m* 


kT 


(- 


2. 

t)3 


(80) 


value  for  A.,  • 
Jk 


'Jk  Jv ik 

T,  . and  T,.  are  the  relaxation  times  of  nuclei  j and  k,  and  the  relaxations 
lj  lk  J 

represented  by  these  times  are  assumed  due  only  to  the  S-character  of 

the  electronic  wave  function  near  the  Fermi  surface;  if  the  electronic 

P-character  contributes  to  the  relaxation  (80)  estimates  too  large  a 

To  obtain  ^80)  Ruder  man  and  Kittel  also  assume  that  the 

electronic  energy  contours  in  k space  are  spherical,  corresponding  to  an 

effective  mass  m*  and  a wave  number  k at  the  Fermi  level. 

m 

. —1  ^ 

For  m -m  and  the  values  of  T ^ in  Table  1,  (80)  predicts  that  h 

SO^for  nearest  neighbors  in  aluminum  and  h A.,  = 140  ^ for  nearest 
63 

Cu  neighbors  in  copper.  These  values  lead  to  additional  contributions 

to  6H  due  to  nearest  neighbors  only  of  about  . 3 gauss  for  both  aluminum 

and  copper.  To  agree  with  experiment  we  have  to  assume  that  Aj^.  is 

about  2 times  its  theoretical  value  in  aluminum,  and  about  3 times  in 

copper.  These  factors  are  rather  large,  but  might  result  from  a small 

effective  mass  or  a complex  band  structure.  An  exchange-type  interaction 

of  this  magnitude  might  still  be  small  enough  to  cause  relatively  little 

exchange  narrowing  of  the  unsaturated  resonance,  as  is  observed.  Un- 

39 

fortunately  for  this  interpretation,  Kambe  has  predicted  theoretically 

sk  **  ^ 

that  in  copper  m = m so  that  an  anomalously  large  A.,  is  not  expected. 

JK 

We  may  summarize  this  section  by  saying  that  the  theory  agrees 
fairly  well  with  experiment  provided  we  assume  that  the  exchange -type 
coupling  constant  A^  is  two  or  three  times  greater  than  its  theoretical 
value,  and  that  the  relaxation  of  neighboring  nuclear  spins  is  coherent 
to  a degree  considerably  greater  than  would  be  expected  from  simple 
considerations.  These  observations  are  consistent  with  the  assumption 
of  an  almost  filled  or  almost  empty  band  with  a small  effective  mass. 
There  appears  to  be  no  independent  support  for  such  an  assumption.  If 
we  assume  instead  a simple  electronic  structure  with  one  electron  per 
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atom  and  m*=m,  the  agreement  between  theory  and  experiment  is  not  so 

4 5 6 

good,  but  is  considerably  improved  over  that  of  previous  * * theories 

of  saturation  as  applied  to  solids.  Finally,  we  should  re-emphasize  the 
fact  that  the  interpretation  of  this  section  is  based  on  a questionable 
calibration,  owing  to  the  disagreement  of  the  proton  calibration  with  the 
other  data . 


V.  Rotary  Saturation 

The  effect  which  we  have  called  rotary  saturation  was  suggested  by 
the  previous  use  of  the  rotating  coordinate  system  representation.  Much 
of  the  theory  above  will  be  useful  in  treating  this  effect  in  solids,  but  it  is 
not  necessary  to  understand  Section  III  in  order  to  understand  rotary 
saturation,  as  it  is  also  a consequence  of  the  Bloch  equations. 

The  effect  is  observed  in  a liquid  obeying  the  Bloch  equations  as 

follows:  the  dispersion  derivative  signal  is  observed  at  resonance  well 

above  saturation,  using  a magnetic  field  modulation  amplitude  which 

is  a sizable  fraction  of  the  rf  amplitude  Hj,  and  a modulation  period 

much  larger  than  T,  and  T7  (or  T , in  the  case  of  a solid).  An  audio 

i.  £.  cg.  4Q 

frequency  magnetic  field  of  frequency  <*>  and  amplitude  H , oriented 

S cl 

in  the  z direction,  is  also  applied  to  the  sample.  When  the  audio  frequency 
approaches  the  frequency 

«a  - YHer  * yH1  (81) 

the  dispersion  signal  is  observed  to  decrease,  and  it  goes  through  a 
minimum  when  the  condition  (81)  is  satisfied.  The  frequency  yHer  is 
the  classical  nutation  frequency  of  the  nuclei  in  the  rf  magnetic  field. 

Figure  9 shows  a rotary  saturation  run  obtained  with  water  heavily 
doped  with  paramagnetic  (Mn++)  ions.  The  dispersion  minimum  frequency 
of  12.75  kc.  corresponds  to  the  proton  resonance  frequency  in  a field  of 
3.0  gauss.  Search  coil  measurements  indicated  that  in  the  run  of  Fig.  9, 

Hj  was  3 gauss  within  the  probable  experimental  error;  these  measurements 
of  were  rather  inaccurate  because  of  the  uncertain  geometric  factors 
involved.  Actually  the  run  of  Fig.  9 was  used  to  calibrate  the  rf  field  for 


r 
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use  in  the  other  runs  reported  in  this  paper,  and  the  rf  field  in  this  run 
was  therefore  assumed  to  be  3.0  gauss.  The  theoretical  justification 
for  this  assumption  will  be  given  below. 


Rotary  saturation  can  be  understood  by  transforming  to  the  rotating 
coordinate  system.  In  Fig.  10a  is  shown  Her  during  a positive  peak  of 
the  14  ~ magnetic  field  modulation.  The  Bloch  equations  predict  that  in 
the  absence  of  the  applied  audio  field  and  in  the  limit  of  large  the 


magnetization  will  be  approximately  in  the  direction  of with  a 


magnitude  M = M cos  © T,T»/T1tI.  At  the  negative  peak  of  the 
Zp  O wl  JLJri 


modulation  the  situation  depicted  by  the  dotted  arrows  in  Fig.  10a 
applies.  The  observed  rf  dispersion  signal  is  proportional  to  M 


zp 


The  action  of  the  audio  field  Ha  in  the  rotating  coordiante  system  can 


be  seen  by  analogy  to  ordinary  saturation  in  the  fixed  coordinate  system. 


Fig.  10b.  If  ®— 90  the  correspondences  H 


er 


H , M — , and  H - 
o op  o a 


2Hj  apply,  and  the  spin  relaxation  processes  are  almost  the  same  in  the 


two  coordinate  systems  (identical  if  T j = T^.  as  is  frequently  the  case). 
In  ordinary  saturation  the  rf  field  2H^  reduces  the  amplitude  of  M ^ to 
some  value  less  than  its  equilibrium  value  M . In  rotary  saturation  the 


audio  field  H reduces  the  amplitude  of  M to  some  value  less  than  its 
a z p 

quasi -static  equilibrium  value  Mq^.  -As  a consequence  the  observed 


dispersion  signal  is  correspondingly  reduced,  as  in  Fig.  9.  Power  is  also 
absorbed  from  the  audio  field,  but  this  absorption  is  normally  too  small 
to  observe  directly. 

Rotary  saturation  can  also  occur  off  resonance,  for  |hq-  y | — Hj, 
although  it  is  most  directly  observed  at  resonance.  As  before,  the  audio 


field  is  most  effective  in  reducing  the  magnetization  M and  the  dis- 
persion signal  when  the  condition  wa  = v^er  obeyed.  The  effectiveness 
of  the  audio  field  is  reduced  by  a factor  sin^  @ since  only  the  square  of 


the  component  of  perpendicular  to  Her  acts  in  the  double  saturation. 


Some  rather  complicated  dispersion  derivative  traces  can  be  obtained  by 
using  various  fixed  values  of«a  greater  than  yH^.  These  signals  all 
appear  to  be  consistent  with  this  simple  picture  . 


Rotary  saturation  can  be  treated  theoretically  in  liquids  by  transforming 


t — r r 


AUDIO  FREQUENCY  (KC.) 


Fig.  9.  Rotary  saturation  run  on  protons  in  water  heavily  doped  with 
manganese  ions.  Hj  was  thought  to  be  3.0  gauss,  Ha  was  about  .6  gauss, 
and  the  magnetic  field  modulation  was  1 gauss  peak-to-peak. 


r 


TR206  -45- 

24 

to  a doubly  rotating  coordinate  system  , the  first  rotation  being  about 
HQ  at  frequency  w,  the  second  about  Her  at  frequency  u>a.  We  ignore 
the  circularly  polarized  component  of  H which  rotates  at  frequency 
2co  in  the  doubly  rotating  system,  and  we  also  ignore  the  time  dependent 

3> 

parts  of  the  relaxation  terms  (involving  and  T£),  which  vary  sinu- 
soidally with  frequency  <oa*  The  resulting  expression  for  is  similar 

to  that  for  as  predicted  by  the  Bloch  equations  for  ordinary  saturation, 
with  the  substitutions  (for  ® — 90°)  M — ^ M , H — »-H  , 2H.— *-H  , M — *-M  , 
oo — ► cvy  ^^2’  an<^  — i *1/2  (Tj1  + T £ * ).  Thus  the  minimum  of  M^p 

and  the  dispersion  signal  should  indeed  occur  when  = yH^,  which 
near  resonance  corresponds  to  = yH^. 

The  time  dependent  parts  of  the  transformed  Bloch  equations  will 
introduce  some  error  in  this  solution  and  are  probably  responsible  for 
the  asymmetry  in  the  run  of  Fig.  9.  They  will  probably  not  introduce 
more  than  a few  per  cent  error  in  the  determination  of  Hj  by  this  method. 

In  Fig.  11  are  shown  some  rotary  saturation  runs  in  aluminum,  for 
various  values  of  In  solids  the  line  width  is  determined  by  the  static 

dipolar  interaction  among  different  spins  rather  than  by  collision  broaden- 
ing as  in  the  case  of  liquids . For  low  rf  fields  the  dipolar  interaction 
displaces  the  minimum  to  a frequency  above  yHj  (indicated  by  arrows  in 
Fig.  11),  but  for  large  the  minimum  very  nearly  coincides  with  the 
predicted  nutation  frequency. 

Rotary  saturation  in  solids  can  be  treated  theoretically  in  the  rotating 

system  in  the  same  way  that  ordinary  saturation  was  treated  by  BPP, 

because  the  analogy  of  Fig.  10  also  applies  to  the  case  of  solids.  The 
29 

effective  energy  absorbed  per  second  by  the  spins  from  the  audio  field 
is  (for  0—90°) 


l/Z^X"  H I 

a a 

where  % "is  the  imaginary  part  of  the  audio  susceptibility  in  the  z,  or 
a 

xp  (approximately)  direction.  Since  the  system  is  assumed  to  be  in  a 

canonical  distribution  of  states  with  respect  to  Jnp  , for  small  H , we 

r er  a 

can  write 


> 

■-  ij 
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=(Uj»(Mcp/Her)f(wa)  (83) 

where  Mop/Her  is  the  rotary  analogue  of  XQ,  and  f(va)  is  a normalized 
shape  function  similar  to  that  defined  by  Broer4l,42.  (83)  is  consistent 
with  a rotary  version  of  the  Kronig-Kramers  relations  and  can  be  derived 
in  the  same  way  as  the  analogous  equation  for  ordinary  paramagnetic 
resonance . 

The  effective  energy  absorbed  per  unit  time  by  the  lattice  from  the 
spin  system  is  (from  Eqs.  37,  47,  and  48) 

t4r  lOir*  -tOir*  o>'Tlp 

SJL.  (84) 

= (M  -M  )H  (1+  6H2/H2J/T, 

' z p op  er  er  lp 

Here<<Je£  ^ = M H (1+  fiH  /H  ),  the  quasi -equilibrium  value  of 

«?6-e  r»* 

In  the  quasi-steady  state  the  thermal  rate  of  decrease  of 
given  by  (84)  must  equal  the  rate  at  which  (effective)  energy  is  absorbed 
from  the  audio  magnetic  field,  given  by  (82).  If  we  assume,  following  BPP, 
that  (83)  holds  under  the  substitution  M — ► M in  arbitrarily  large  H , 

— Z p 3l 


we  get 


M /M  n = [1 +0/2)H2«2f(v  )T.  H"2(l  + 6H2/H2  l"1]"1  (85) 

zp  op  L ' ’ a a ' a'  lp  er'  er  J 

The  observed  rf  susceptibility  is  thus  reduced  oy  a factor  equal  to  the  right- 

hand  side  of  (85). 

Actually,-' the  assumption  that  M^p  can  substituted  for  Mq  in  (83)» 

for  arbitrary  H , is  not  justified.  This  is  the  same  erroneous  assmption 
a 5 

made  in  the  BPP  theory  of  ordinary  saturation,  as  discussed  in  Section 
111(a).  Equation  (85)  is  expected  to  hold  only  for  small  H , where 

Mzp/H>p  -1- 

The  shape  function  f(v  ) can  be  determined  experimentally  by  making 

3l 

a series  of  runs  at  constant  H,  and  different  H and  w , and  using  (85)  to 

determine  f(v  )T , (1+  fiH  /!*_)“  ^.  Data  obtained  in  this  way  for  aluminum 
& ip  cr 

and  copper  are  shown  in  "Figs.  12  and  13.  In  drawing  the  experimental 
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curve  (solid  line)  through  these  data  we  have  given  the  greatest  weight 
to  those  points  obtained  at  low  audio  power  levels,  since  (85)  is  expected 
to  be  more  accurate  at  these  levels. 

Integration  of  the  solid  curves  in  Figs.  12  and  13  yields  values  of 
2 2 _ 1 

T,  (1+6H  /H.)  of  5.6  milliseconds  for  copper  and  2.2  milliseconds 

l p i 

for  aluminum.  These  times  differ  from  values  predicted  from  the 
previously  estimated  and  6H  by  a factor  of  about  two.  The  value 


for  copper  is  too  large  and  that  for  aluminum  is  too  small, 
for  this  discrepancy  is  not  known. 


The  reason 


We  have  not  attempted  to  analyze  the  shape  function  f(v  ) in  detail, 
but  we  can  make  a few  general  remarks.  A considerable  number  of 

theoretical  predictions  conerning  f(v  ) can  be  made  using  the  methods 

25  42  a -if 

of  Van  Vleck  and  Wright,  but  the  Hamiltonian  J°tT  has  fewer  secular 

terms  than  the  spin  Hamiltonian  in  the  fixed  coordinate  system,  and 

® 42 

therefore  the  theoretical  work  of  Wright  is  not  directly  applicable  here. 
As  expected,  f(v)a  has  definite  peaks  at  zero  frequency  and  at  yH^.  The 
significance,  if  any,  of  the  small  departures  of  the  resonant  peaks  from 
the  value  yH^  is  not  known.  These  departures  may  be  due  in  part  to 
experimental  error. 

In  the  limit  of  large  H , f(v  ) is  expected  to  be  a symmetrical  curve 

cr  e 

centered  about  the  resonance  frequency  yH.,  and  the  moments  of  f(v  ) 

X cL 

in  this  case  have  been  calculated  by  Van  Vleck.  It  is  noteworthy  that  if 

A..  = 0 the  secular  part  of  (first  three  terms  of  Eq.  20)  contains  a 

j k e p 

spin  interaction  identical  in  form  and  exactly  half  as  large  as  the  secular 

part  of  the  Hamiltonian  in  the  fixed  coordinate  system.  This  implies  that 

g(v),  the  rf  line  shape,  should  be  exactly  similar  and  twice  as  broad  as 

f(v  ).  The  expected  shapes  f(v  )~g(2v)are  plotted  in  Figs . 12  and  13; 

here  g(v)  was  obtained  experimentally.  The  observed  f(v  ) is  considerably 

narrower  than  the  predicted  shape,  which  may  be  due  in  part  to  the  fact 

that  Hj  is  not  very  large  and  is  in  fact  comparable  to  the  predicted  line 

width.  Another  possible  reason  for  the  narrow  observed  f(v  ) may  be 
~ a 


the  terms  A.,I.'I,  , which  are  expected  to  produce  exchange  narrowing  of 
f(v  ).  It  is  rather  surprising  that  they  do  not  also  produce  appreciable 

cL 

exchange  narrowing  of  the  rf  resonance  line. 
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It  is  amusing  to  note  that  if  l-3cos2  ® = 0 and  Hj  is  large,  the  line 
width  should  be  zero,  except  for  a small  spin-lattice  relaxation  broaden- 
ing. In  the  case  of  the  data  of  Figs.  12  and  13  the  modulation  was  not 
large  enough  to  induce  such  narrowing. 

In  conclusion  we  may  say  that  besides  providing  an  accurate  calibration 
of  the  rf  field,  rotary  saturation  is  a relatively  simple  way  to  study  the 
properties  of  spin  systems  in  small  magnetic  fields.  This  statement  is 
based  on  the  fact  that  the  time  dependent  part  of  Z*r  really  can  be 
neglected  to  a very  good  approximation,  so  that  rotary  saturation  is 
really  closely  analogous  to  ordinary  saturation. 


VI.  Concluding  Remarks 

The  theory  and  experiment  in  this  paper  demonstrate  the  usefulness 
of  the  rotating  coordinate  representation  in  resonance  problems.  The 
rotating  coordinate  representation  should  be  useful  in  treating  other  types 
of  relaxation  and  spin  interaction  than  those  considered  here.  For  example, 
the  rotary  saturation  experiment  indicates  that  the  perturbations  re- 
sponsible for  transverse  (T 2 or  T 2&)  relaxation  of  nuclei  are  those  of 
frequency  yHj,  rather  than  zero  frequency. 

All  the  theory  in  this  paper  can,  of  course,  be  applied  to  paramagnetic 
resonance  under  suitable  conditions.  In  the  case  of  ferromagnetic  resonance 
the  approximations  made  above  are  presumably  invalid,  but  some  progress 
might  be  made  along  similar  lines. 

An  as  yet  unexploited  consequence  of  this  work  is  the  feasibility  of 
measuring  spin  temperatures  and  relaxation  times  and  obtaining  nega- 
tive temperatures  in  solids  by  the  method  of  adiabatic  fast  passage.  The 
experimental,  technique  would  be  similar  to  that  employed  by  Drain43  and 
others  to  measure  relaxation  times  in  liquids. 

The  implications  of  this  research  concerning  the  interpretation  of 
previous  work  are  rather  unimportant,  except  for  the  probable  errors  in 
previous  fast  modulation  saturation  measurements  of  Tj,  as  discussed  in 
Section  III.  The  optimum  signal  to  noise  ratio  predicted  by  Bloch  and 

3??  iS  t0°  SmaU  by  a factor  of  T2/Tl,  provided  the  dispersion  mode 
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is  observed  under  conditions  of  slow  passage.  This  gain  is  usually 
offset  by  the  additional  apparatus  noise  at  high  rf  power  levels,  and 
by  the  fact  that  Tj  is  frequently  so  long  that  slow  passage  is  impractical. 
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APPENDIX  A 


The  method  used  to  calculate  the  ratio  of  the  effective  external  energy 

M H to  internal  spin-spin  effective  energy  is  similar  to  that  .used  by 
zp  er 

Van  Vleck  to  calculate  the  specific  heat  of  a spin  system.  However, 

the  form  of  the  Hamiltonian  is  slightly  different  and  the  calculation  is 

21 

simplified  by  the  use  of  the  density  matrix  formalism. 

The  assumption  of  a canonical  distribution  of  states  with  respect  to 
the  transformed  spin  Hamiltonian  is  equivalent  to  the  assumption  that  the 
state  of  the  solid  is  decribed  by  the  density  matrix  (in  the  r-rotating 
system) 


P^Cexp  (-#er/kT*  ->^/kT). 


(86) 


Here  C is  a normalizing  factor  making  Trp  = 1,  T*  is  a constant  analo- 
gous to  an  effective  spin  temperature  in  the  rotating  system, is  the 
lattice  (non-spin)  Hamiltonian,  and  T is  the  lattice  temperature.  TV  O 
denotes  the  trace  (diagonal  sum)  of  the  operator  O.  The  spin-lattice 
interaction  J 

is  neglected  since  it  is  assumed  to  be  small.  The  first  term  of  the  ex- 
ponential (86)  represents  the  assumption  that  the  spin  system  is  in  a 

canonical  distribution  of  states  with  respect  to  , corresponding  to 

6 r 

the  spin  effective  temperature  T . The  second  term  represents  the 
assumption  that  the  lattice  (electrons,  lattice  vibrations,  etc.)  is 
relatively  unaffected  by  the  rf  field  and  can  be  described  by  the  temper- 
ature T . 


(or  transformed  to  the  rotating  coordinate  system) 


In  practice  we  can  assume  that  T 


* 


is  large  enough  for  us  to  write 


P — ( 1 ->4r/kT*)pL, 


(87) 


where 


PL=  C exp  ( -J^/kT). 


(88) 


The  expectation  value  averaged  over  the  canonical  ensemble 


is 
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«J4  »^Tr  J4  p = (TtU  p - TrK2  p /kT*  ) . 
er  err  errL  err  L ' 


(89) 


This  expression  is  easily  evaluated  using  a representation  in  which 

the  operators  I.  a.ndJ=fT  are  diagonal.  In  this  case  we  can  write 
- 1 n ^ 2 if 

C = (2I  + 1)  Tr ^exp( -/=tj^/kT)  where  N is  the  total  number  of  spins 
and  Tr^O  denotes  the  diagonal  € am  of  the  operator  O over  all  the 
eigenvalues  <5^  of  keeping  the  quantum  numbers  m^ . (eigenvalues 


of  Ijz)  fixed.  In  the  same  representation  we  can  write 

=<Tri^r»Tri.Pi J ■ 


(90) 


2 2 

where  TrTHrer  denotes  the  diagonal  sum  ofj£=£er  over  the  quantum  numbers 

rrijj,  keeping  the  quantum  number  fixed.  This  relation  follows  directly 

from  the  diagonality  of  with  respect  to  the<5^  and  of«^  with  respect  to 

the  mT . . 

Ij 

Using  the  easily  verified  relations 

TrI  W'  >»I(I+l)(2I  + 1,N6j  k6  ^ , 

* 1/9  l2<»l>2UI+»t\k<V  v6m,n6W  (92) 
where  v ,p.  ,Y^  = x,y,z  and  j ^ m,  k ^ n, 

we  can  evaluate  TrJ=(?rp  ^ an<i  similarly  show  that  TrJr  rp^  is  zero.  The 
result  is 


(91) 


-kl*«>£er»=  1/3  Ng2p2H2r  1(1+1)  + 1/9  I2(I+1)2 


Z15^*2'3  Bj2k'- 


(93) 


The  first  term  o£<&Ja&  y>  is  the  external  energy-M  • H , as  can  be 
verified  by  directly  calculating  « Mr»=  «g{3 using  (87).  The 
second  term  is  the  internal  spin- spin  effectivi  energy.  Equations  (30) 
and  (31)  follow  directly  from  (93)  on  comparison  with  the  Van  Vleck2^ 
expression  for 
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APPENDIX  B 


The  equation  of  motion  of<tJe£er»  is 


— «Wer»  ' «!  Wr'^erl^  + « 


8*lr 


(94) 


dt 


at 


where  J£r  is  the  total  Hamiltonian  + -^g  trans*orme(*  to  t^ie 

r -system; 

*r=-*L+-Wt)+>tr-  <95> 

Since  -^r  commutes  with  itself  and^^,  (94)  becomes 


dt 


<<7H  >>  = 4:9^  >>— M • H • 

"'rcer'  er  wr  «ver 


(96) 


The  abbreviation  90  is  used  for  -ih  30  t6at<OO^is  the  time 

rate  of  change  of  the  canonical  average  expectation  value  «0»of  0,  due 
to  the  (time  dependent)  spin-lattice  interaction  ^ . In  the  treatment 
of  fast  or  intermediate  passage  the  first  and  last  terms  of  (96)  are  re- 
tained ( Eq.  (65),  Section  111(e)  ),  while  under  steady-state  conditions  they 
are  by  definition  zero,  and  the  second  term  iB  therefore  also  zero. 

«8jfr»  is  also  referred  to  in  Section  III  as  [8/8t]g j5<^er^”  Assuming 

that  the  spin-lattice  interaction and  therefore  the  transformed 

interaction , is  linear  in  the  operators  I.  , I.  , I.  we  get 
SLr  r jx’  jy*  jz 

2_  2 


«aXer»  = Tr (01.)  p+  g P ^[Aj2k  Tr(81.)-Ikp 


(97) 


+ BjkTr(a  V^- 

In  (97)  P is  understood  to  be  given  by  (87). 

To  evaluate  (97)  we  use  the  assumption  (34),  which  in  slightly  more 
general  form  can  be  written 

P'/T,_»  (98a) 


Tr  (81.  )p'  = - Tr  I. 

jx,r  jx 


2e 


Tr(8Ijy)P'  = - Tr  I.y  P*/T2e 


(98b) 
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TrOI.jp*  = -(TrIjzpl-Io)/T1. 


(98c) 


Here  we  have  temporarily  abandoned  the  assumption  Tj  = T2g,  and  p1  is 
a density  matrix  describing  an  arbitrary  state  of  the  spin  system.  Thus 
p’  = pj  pj^,  where  pj  is  an  arbitrary  function  of  the  spin  operators  for  which 
Trj  Pj  = (21+1)^.  Use  of  the  various  operators  1,  1+1^,  1+1^1^,  1+1^ 

1^  (with  ^ ,1}  , £=  x,  y,  z and  j,  k and%  ,7)  , £ not  necessarily  identical 
or  non-identical)  for  pj  yields 


Tr<aijz»L 

VT1  ' 

(99) 

Tr«VVi.  = 

- 1/3  1(1+  1)/T2e  . 

(100) 

Tr<8Ijz)Ij.PL  ' 

- 1/3I(I+1)/Tj  . 

(10D 

Tr(9I.  )I.  l£  p = 
jz'  jz  kv  L 

- 1/9  I2(I  + 1)2/TX  , 

(102) 

Tr(8I.  )I.  I.2  pT  = 
jv'  jv  k|irLi 

- 1/9  I2(I+  1)2/T2e  , 

(103) 

s*  v = x,  y.  All  other  traces  occurring  in  (97) 

are 

zero.  Using  these  expressions  we  get 


2 02, 


«»H,8  * Ng«H0-  — U0/T1  - N*  f[!S!1»tl>[(Ho-  — )2/rltH^/T2e  ] 


3kT 


!Z(1-H)2XI  +-g?—  ] 1 ■ (104) 

jfk  I Jk  ^ ^ ?T1  9TX  J 

«Mr»  =Trgp^Lp  is  easily  shown  to  be  a vector  in  the  direction  of 
> 4 J 


H or  z , of  magnitude 
^er  p ° 


M = I(I+l)gZpZNHer 
Z P 3kT* 


(103) 


The  static  equilibrium  magnetization  Mq  is  given  by  the  same  formula  with 
Mq,  T->T,  andHer— >-Ho>  Therefore  (104)  can  be  written 


«8>tr»=  . (M  - Mop>(l+  6H2/H^r)  Hgr /T  . 


(106) 
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where 


2/„2 


><1+  6H  /Her^l 


M cos  ® 
o 


M cos  ®(Tr  ^os2  0+ Tl1  sin2  0)  1 
o ' 1 2e 


Tj(^a  6H7H*r) 


UAj^0,  a^-^  + glj-Kf^cos2®  + T-‘  .i 


. 1 2 “ 1 
2e  8in  ®> 


(107) 


(108) 


If  T^e  = Tj,  and  A^  is  not  necessarily  zero,  Eqs.  (38)  and  (63)  of 
Section  III  follow  immediately. 

In  this  calculation  and  cor  tain  only  the  time  dependent  part 

of  the  spin-lattice  interaction  which  is  responsible  for  the  relaxation,  while 

a-ndjzf  contain  time  independent  parts  which  result  from  the  spin-lattic* 
o er 

interaction  (i.e.,  chemical  shift,  Knight  shift,  and  electron  coupled  spin- 
26  27 

spin  interactions  * ).  As  discussed  in  Section  III(c),  this  procedure  is 

hard  to  justify  rigorously  but  seems  reasonable  physically. 


APPENDIX  C 


We  assume  that  the  system  is  decribed  by  the  density  matrix 


p = C exp 


r%rl/kT%  [ gpHo  s iyt*  hjm] 


(109) 


wherej^grj  contains  all  the  time-independent  terms  ofJ^r  (Eq.  (53)  ) ex- 
cept the  last.  The  reason  for  including  aj.1  the  spin-spin  terms  of  is 

discussed  in  Section  111(d). 

The  ratio  of  external  unprimed  spin  energy  to  internal  spin-spin  energy 
is  obtained  by  evaluating^J=^r  Proceeding  as  in  Appendix  A,  (55) 


} 

- —-eZj 
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easily  follows . 

The  equation  of  motion  o f^>^£  is  given  by  (96)  with  the  substitutions 

J=^er and  die  density  matrix  (109)  for  that  of  (86).  Evaluation  of 

the  various  terms  is  similar  to  the  procedure  of  Appendix  B.  Equations 
(99)  to  (103)  hold  for  the  substitutions  j— *-j'  and/or  k— >k'  under  suitable 
additional  substitution  ofJ^j,  Tj,  and  T^g.  (48)  and  (56)  then  follow  easily 
from  the  assumption  and  T'j  = ^2e* 
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